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ABSTRACT 

'  a  ’ 

The  purpose  of  this  research-' was  to  examine  robustness, 
mean  relative  tightness,  coverage,  and  power  of  selected 
unmodified  and  modified  Cox  and  Snell  and  Stringer  error 
limit  bounds.  The  simulation  was  performed  by  repetitive 
sampling  from  an  accounting  population  with  various  known 
book  value  and  error  distributions.  Additional 
modifications  to  the  selected  modified  Cox  and  Snell  bounds 
was  done  by  incrementally  loosening  the  bounds  by  5,  10,  and 
20  percent  in  a  search  for  bounds  with  better  performance 
characteristics . 

There  were  several  conclusions  that  could  be  made  from 
this  research.  The  modified  Cox  and  Snell  bounds  can 
achieve  high  coverages  for  accounting  populations  with  low 
error  amount  intensity  (EAI )  with  significant  increases  in 
mean  relative  tightness  over  unmodified  Cox  and  Snell 
bounds.  The  Stringer  and  Cox  and  Snell  bounds  can  still 
achieve  high  coverages  with  significant  improvements  in  mean 
relative  tightness  when  the  nominal  confidence  level  is 
lowered  from  90%  to  85%.  Only  minor  changes  in  prior 
probability  settings  materially  affect  the  performance  of 
Cox  and  Snell  bounds.  And,  in  accounting  populations  with 
low  EAI,  the  selected  Cox  and  Snell  bounds  are  conservative. 


ERROR  PREDICTIONS  IN  ACCOUNTING  POPULATIONS 


I .  Introduction 

General  Problem 

Auditors  are  required  to  make  statements  about  the 
accuracy  of  financial  accounts  by  reviewing  internal 
controls  and  empirical  sampling.  Through  sampling  and 
inferential  statistics,  auditors  estimate  as  accurately  as 
possible  tne  maximum  amount  an  account  is  likely  to  be  in 
error.  Their  estimate  is  known  as  the  upper  error  limit. 

The  upper  error  limit  is  then  compared  to  the  material 
error,  which  is  the  limit  of  precision  necessary  for  the 
account.  Results  of  the  comparison  enable  the  auditor  to 
state  an  opinion  about  the  reasonable  accuracy  of  the 
account  (8:4,5) . 

The  estimation  of  upper  error  limits  depends  on 
our  ability  to  construct  models  incorporating  two 
aspects  of  accounting:  the  error  rate  and  the 
distribution  of  error  sizes  for  those  accounts  in 
error  £3:1 80  J . 

Traditionally,  to  reduce  risk,  auditors  have  chosen  upper 
error  limit  bounds  which  result  in  very  broad  confidence 
intervals,  possibly  far  in  excess  of  the  actual  error.  This 
conservative  approach  is  due  to  the  fact  that  little  is 
known  about  the  characteristics  of  accounting  populations 
(7:270).  Conservatism  leads  to  possible  erroneous  rejection 
of  the  hypothesis  that  the  book  value  of  the  audited 


population  is  reasonably  accurate.  Therefore,  it  is 
necessary  to  choose  sampling  plans  and  error  limit  bounds 
which  are  appropriate  to  the  auditing  situation. 

Sampling  should  be  done  in  such  a  way  that  a  maximum 
amount  of  information  can  be  obtained  at  a  minimum  cost. 

When  the  auditor  samples  accounting  populations, 
he  encounters  two  key  sets  of  problems.  One  is 
connected  with  low  error  rates  in  the  accounting 
populations.  The  otner  pertains  to  the  effectiveness 
of  sampling  procedures  designed  for  either  low  or  high 
error  rates  when  the  actual  error  rates  are  not  of  the 
anticipated  magnitude  £  1 0  :  4  J . 

Low  error  rates  pose  significant  problems  for  auditors  who 
use  classical  statistical  estimators.  Dollar-unit  sampling 
(DUS)  was  "specifically  developed  for  use  in  the  auditing 
environment  where  low  error  rates  in  account  balances  are 
often  encountered"  (3:179).  Bayesian  procedures  using  prior 
subjective  information  about  accounting  populations  have 
been  used  DUS  to  produce  tighter  upper  error  limits  than 
those  obtained  with  the  DUS  method  and  classical  estimators 
(3:179).  The  Bayesian  approach,  though,  has  its  own 
problems:  (1)  difficulty  in  assessing  a  prior  probability 

judgment,  and  (2)  measurement  of  utility  in  those  decisions 
where  the  sums  involved  are  large  (1:115). 


Definition  of  Key  Terms 


Confidence  Interval: 


the  range  around  a  calculated 


sample  mean  in  which  the  population  mean  is  expected  to  lie 
at  tne  given  confidence  level. 


Nominal  Confidence  Level:  the  expected  frequency  of 
correct  intervals  or  bounds  in  repeated  samples. 

Conservatism:  larger  upper  error  limits  calculated 
than  the  nominal  confidence  level. 

Error  Taint:  the  error  amount  found  in  a  sample  unit 
divided  by  the  true  book  value. 

Rel  Tightness:  the  magnitude  of  an  upper  error  limit 
relative  to  the  true  error  in  the  population. 

Coverage:  in  repetitive  sampling,  the  percentage  of 

calculated  upper  error  limits  which  are  equal  to  or  greater 
than  the  true  dollar  error  in  the  population. 


Backg  round 

Upper  Error  Limits.  When  auditors  first  began  using 
sampling  techniques  to  estimate  errors  in  accounting 
populations,  they  used  classical  sampling  estimators:  ratio 
and  difference  estimators. 


Stringer  (1963)  cautioned  that  use  of  these 
estimators  entails  potential  dangers  because  their 
estimated  standard  errors  will  equal  zero  when  no 
errors  are  found  in  the  sample.  With  classical  large- 
sample  confidence  limits  this  would  imply  that  the 
point  estimate  of  the  total  audit  amount  or  the  total 
error  amount  is  perfect,  which,  obviously  is  not  a 
reasonable  conclusion  in  an  audit  situation  employing 
a  moderately  large  sample  size  19:77], 


To  overcome  the  observed  shortcomings  with  ratio  and 
difference  estimators.  Stringer  developed  his  own  method  of 
calculating  error  limits  (called  the  Stringer  bound) . 
Although  the  Stringer  bound  has  been  used  extensively,  it  is 
not  without  its  own  problems.  First  of  all,  no  general 


theory  has  been  published  which  supports  the  Stringer  bound. 
This  lack  of  theory  occurs  because  the  Stringer  bound 
involves  a  combination  of  attributes  and  variables 
principles  which  makes  theoretical  analysis  difficult. 
Second,  the  Stringer  bound  itself  is  very  conservative, 
resulting  in  upper  error  limits  far  greater  than  the  true 
error  in  populations  (9:78). 

In  1973,  Anderson  and  Teitlebaum  also  identified 
potential  dangers  in  using  ratio  and  difference  estimators 
for  auditing  when  the  error  rate  is  small.  They  proposed 
using  dollar-unit  sampling  together  with  the  Stringer  bound 
for  estimating  the  total  overstatement  error  in  accounting 
populations.  Dollar-unit  sampling,  which  considers  the 
individual  dollar  as  the  sample  unit,  requires  no  further 
stratification  by  book  amount  since  all  sampling  units  are 
exactly  the  same  size  in  terms  of  book  value  without  the 
need  for  stratification  [8:78, 79 J.  "Dollar-unit  sampling 
can  be  thought  of  as  employing  the  ultimate  in 
stratification  by  book  amount"  [9:78J. 

In  order  to  obtain  more  empirical  evidence  about  the 
behavior  of  classical  estimators,  Neter  and  Loebbecke  in 
1975  conducted  simulation  studies  with  four  actual 
accounting  populations.  They  examined  the  performance  of 
the  raean-per-unit ,  ratio,  and  difference  estimators  each 
with  simple,  stratified  random,  and  dollar-unit  sampling. 
They  found  that,  for  the  estimators  under  study,  performance 


varied  significantly  depending  upon  the  error  sizes  and  the 
distribution  of  errors  within  the  populations.  When  the 
population  error  rate  was  small  to  moderate  or  the 
distribution  of  errors  was  highly  skewed,  the  actual 
proportions  of  correct  confidence  intervals  was  far  below 
the  nominal  level.  Even  under  the  most  favorable 
conditions,  Neter  and  Loebbecke  found  that  there  is  the 
possibility  that  the  nominal  confidence  coefficient  based  on 
the  normal  distribution  differs  somewhat  from  the  actual 
confidence  level  110:134,135]. 

In  1977  Neter,  Leitcn,  and  Fienberg  proposed  a 
different  method  of  obtaining  a  oound  for  total  population 
error.  Basing  their  bound  on  the  multinomial  probability 
distribution,  they  named  it  the  multinomial  bound  for 
overstatement  errors.  The  multinomial  bound  has  been  found 
to  be  effective  because : 

( 1 )  it  incorporates  information  about  the  magnitude  of 
the  sample  errors;  (2)  it  has  distributional 
properties  which  are  fully  known;  and  (3)  it  does  not 
require  any  assumptions  about  the  distribution  of  the 
population  or  of  the  errors  in  the  population  [5:3]. 

When  used  with  dollar-unit  sampling,  the  multinomial  bound 

has  been  shown  to  be  significantly  tighter  than  the 

Stringer  bound,  with  actual  confidence  levels  nearer  to 

nominal  levels  (5:3). 

While  the  use  of  monetary-unit  sampling  with 
either  the  Stringer  bound  or  the  multinomial  bound 
avoids  the  problem  of  actual  confidence  levels  far 
below  the  nominal  level  when  there  are  few  errors  in 
an  accounting  population,  neither  the  Stringer  bound 
nor  the  multinomial  bound  permit  a  direct  netting  of 


over-  and  understatement  errors.  Further,  neither  of 
these  two  bounds  directly  takes  into  account  prior 
information  which  the  auditor  may  have  about  the 
population  error  rate  and/or  the  magnitude  of  the 
taintings  in  the  population. 

In  contrast,  a  Bayesian  approach  to  constructing 
bounds  for  the  total  error  amount  in  an  accounting 
population  does  enable  the  auditor  to  incorporate 
prior  information  into  the  estimation  procedure, 
though  there  may  be  difficulties  in  making  the  needed 
assessments.  The  Bayesian  approach  also  can  provide 
bounds  not  only  for  the  total  overstatement  and/or 
understatement  error  amounts  in  an  accounting 
population  separately,  but  for  the  net  total  error 
amount  as  well  [5:3]. 

A  great  deal  of  support  has  been  generated  to 
integrate  prior  subjective  knowledge  of  accounting 
populations  into  statistical  calculations  from 
statisticians  as  well  as  auditors. 

Thus  perhaps  the  strongest  argument  favoring  the 
formal  Bayesian  decision  process  is  that  this 
intuition  [prior  subjective  knowledge]  is  incorporated 
formally  into  the  decision  process.  Moreover,  by 
quantifying  the  intuition  in  the  model,  the  Bayesian 
approach  tends  to  prevent  haphazard  use  of  this 
experience  [1:111]. 

My  own  view  on  the  problem  of  estimating  the 
money  value  of  error  is  that  a  Bayesian  approach  is 
the  only  satisfactory  one  from  a  statistical 
viewpoint.  The  auditor's  judgement  on  many  issues  is 
highly  subjective  and  thus  there  is  no  reason  why 
subjective  ideas  should  not  be  introduced  into  the 
interpretation  of  his  test  results  [12:279], 

Cox  and  Snell  (1979)  proposed  the  use  of  "infinite 

population  theory  as  a  theoretical  foundation  for  the 

development  of  a  Bayesian  upper  bound  for  the  total 

overstatement  amount  in  the  population"  (5:499).  A  1984 

study  done  by  Godfrey  and  Neter  produced  the  following 

conclusions  concerning  the  Cox  and  Snell  bounds:  (1)  the 


bounds  are  substantially  affected  oy  prior  parameter 
assessments,  (2)  they  are  much  tighter  than  the  Stringer 
bound  and  have  high  coverage  when  the  actual  population 
contains  low  error  rates,  and  (3)  they  are  moderately 
tighter  than  the  Stringer  pound  for  populations  whose  error 
amount  is  large  (5:520). 

Godfrey  and  Andrews,  in  1982,  presented  an  alternative 
to  classical  and  other  Bayesian  procedures.  Their  method, 
the  finite  Bayesian  procedure,  correctly  assumes  a  finite 
population  (which  other  Bayesian  methods  do  not)  and 
allows  the  integration  of  the  auditor's  prior  information 
through  a  prior  distribution.  They  found  that  their  finite 
Bayesian  model  more  closely  approaches  the  actual 
population  than  other  procedures  and  usually  requires 
smaller  sample  sizes.  As  with  other  Bayesian  models,  the 
difficulty  in  the  finite  Bayesian  procedure  lies  in 
determining  a  reasonably  accurate  prior  assessment  (4:304- 
315)  . 

Several  other  Bayesian  approaches  have  been  proposed 
by  statisticians  and  auditors.  McCray,  in  1982,  assumed  a 
discrete  prior  distribution  for  the  total  error  amount  and 
a  multinomial  distribution  for  the  monetary  unit  taints. 
Garstka  proposed  a  compound  Poisson  model  which  uses  a 
discrete  prior  distribution  for  the  mean  error  taint  which 
is  used  to  calculate  an  upper  error  bound  for  the  total 
error  amount.  Vanacek's  approach  uses  a  beta  prior 


7 


distribution  for  the  population  error  rate,  then  obtains 
Bayesian  posterior  bounds  for  the  population  error  rate 
(5:499) . 

In  1984,  Dworin  and  Grimlund  published  a  paper 
proposing  a  new  method  of  statistical  analysis  for  dollar- 
unit  sampling.  Their  method  was  designed  to  overcome  the 
conservative  nature  of  previously  applied  upper  error 
bounds.  Their  new  bound,  called  the  moment  bound,  holds 
considerable  promise  for  inventory  accounts.  Results  from 
their  study  indicate  that,  when  the  material  limit  of  the 
population  is  less  than  four  percent  of  book  value,  the 
moment  bound  will  lead  to  the  acceptance  of  more  accounts 
than  the  multinomial  bound.  For  populations  where  both 
understatement  and  overstatement  errors  are  present,  the 
moment  bound  provided  significantly  tighter  bounds 
(2:219,23b) . 

In  1985,  Phillips  examined  modified  and  unmodified 
Stringer  and  Cox  and  Snell  bounds  for  robustness  using  96 
accounting  populations.  He  found  that  none  of  the  modified 
Cox  and  Snell  bounds  were  robust  and  that  both  the  modified 
and  unmodified  Stringer  bounds  were  robust  at  the  95  and  85 
percent  levels  of  confidence.  The  modified  Stringer  bound 
produced  upper  error  limits  which  were  approximately  19 
percent  tighter  relative  to  the  unmodified  Stringer  bound 
[11:105,106].  Phillips  also  found  that  decreasing, 
uniform,  and  increasing  clustering  of  errors  among  line 
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items  did  not  have  a  significant  impact  upon  bound 
performance . 

In  a  study  done  by  Helton  in  1985,  modified  and 
unmodified  DUS-cell  and  Stringer  bounds  were  examined  for 
robustness  at  varied  confidence  levels.  For  the 
populations  under  study,  he  found  that  the  only  robust 
bounds  were  the  unmodified  Stringer  bounds.  However,  all 
coverage  failures  occurred  when  a  50  percent  line  item 
error  rate  was  introduced  into  the  study  population.  Also, 
under  certain  error  conditions,  the  unmodified  Stringer  and 
DUS-cell  bounds  produced  tighter  bounds  relative  to  the 
respective  modified  bounds  [6:4-1]. 

Characteristics  of  Accounting  Populations.  Auditors 
require  knowledge  of  the  frequency,  magnitude, 
distribution,  and  the  possible  causes  of  errors  in 
accounting  populations  to  effectively  and  efficiently 
provide  statements  about  the  accuracy  of  the  account 
(7:271).  Prior  to  a  study  done  by  Johnson,  Leitcn,  and 
Neter  in  1981,  little  had  been  published  about  accounting 
population  characteristics.  Their  study  concentrated  on 
accounts  receivable  and  inventory  audits  and  analyzed  the 
magnitude  of  error  rates,  the  susceptibility  of  individual 
items  to  error  within  the  account,  and  the  degree  of 
overstatement  and  understatement  errors  (7:273,274). 
Highlights  from  the  summary  of  their  results  are: 

1.  There  is  great  variability  in  the  error  rates  for 

both  types  of  audits,  with  the  inventory  audits 


tending  to  be  substantially  higher  than  those  for 
accounts  receivable. 

2.  There  is  some  evidence  from  the  81  audits  which 
suggests  that  the  error  rates  in  both  types  of  audits 
may  be  higher  for  the  larger  accounts  and  for  accounts 
with  larger  line  items  than  for  other  accounts. 

3.  Most  errors  in  receivables  audits  are 
overstatement  errors,  wnile  inventory  audits 
overstatement  and  understatement  are  more  balanced  in 
number . 

4.  The  distribution  of  taintings  are  variable  for 
both  types  of  audits,  especially  those  for  inventory 
audits,  and  depart  substantially  from  a  normal 
distribution.  Quite  a  few  of  the  tainting 
distributions  for  inventories  are  negatively  skewed 
[7:291]. 

These  results  verify  that  great  care  must  be  taken  to 
choose  the  appropriate  methods  for  sampling,  estimating 
population  means,  and  calculating  upper  error  limits. 
Assuming  a  normal  distribution  is  not  appropriate  for 
accounting  populations.  Results  from  the  study  indicate 
that  dollar-unit  sampling  is  the  most  beneficial  method,  in 
terms  of  gathering  accurate  information,  of  sampling 
accounting  populations. 

Conclusion.  Auditors  are  faced  with  significant 
problems  in  choosing  appropriate  sampling  techniques  and 
upper  error  limits.  Much  research  has  been  done  trying  to 
deal  with  these  problems,  but  no  clear  solution  currently 
exists.  Further  research  in  examining  current  methods  and 
the  development  of  new  methods  is  needed.  More  efficient 
and  effective  auditing  procedures  will  assist  auditors  in 
making  more  accurate  and  reliable  decisions. 


SPECIFIC  PROBLEM 


Auditors  must  deal  with  two  basic  problems  with 
accounting  populations.  First,  accounting  populations  have 
relatively  few  errors.  Sampling  the  populations  may  result 
in  finding  very  few  errors,  which  poses  significant  dilemmas 
when  dealing  with  classical  statistical  estimators  such  as 
the  mean-per-unit ,  ratio,  or  difference  estimators.  Using 
the  sample  results  to  construct  confidence  intervals  for  the 
population  error  may  be  unacceptable.  The  actual  confidence 
intervals  provided  by  classical  estimators  may  be  far  less 
than  their  corresponding  nominal  levels.  This  leaves  the 
auditor  with  no  statistical  foundation  upon  which  to  base 
his  audit  opinion.  The  second  major  problem  facing  auditors 
is  determining  an  acceptable  method  of  quantifying  the 
judgment  an  auditor  wants  to  make  about  an  account. 

Classical  statistics  provide  no  way  for  auditors  to 
explicitly  incorporate  subjective  opinions  into  the  final 
judgment.  Bayesian  statistical  models  provide  a  method  by 
which  auditors  may  introduce  their  professional  judgments  of 
particular  accounting  populations  formally  into  upper  error 
limit  calculations  and  also  help  deal  with  the  dilemma  of 
low  error  rates.  Proposed  modifications  to  Cox  and  Snell 
Bayesian  models  have  yet  to  be  fully  tested  and  formally 


Research  Objective  1 ♦  The  first  objective  of  this 
research  study  is  to  examine  robustness  at  the  90% 
confidence  level  for  selected  modified  Cox  and  Snell  upper 
error  limit  bounds  using  study  populations  with  known  mean 
book  values,  error  rates,  and  error  magnitude  distributions. 
For  the  particular  model  to  be  robust,  it  should,  with 
repetitive  sampling,  be  correct  at  (1  -o»)  x  100  percent  of 
the  time  for  each  study  population. 

Research  Objective  2.  The  second  research  objective  is 
to  examine  the  bound  performance  using  measurements  of 
coverage  and  rel  tightness.  Coverage  is  the  proportion  of 
the  time,  using  repetitive  sampling,  an  upper  error  limit 
equals  or  exceeds  the  true  error  in  each  study  population. 
Rel  tightness  is  a  ratio  of  one  conservative  bound  to 
another  less  conservative  bound.  For  the  purpose  of  this 
research,  all  modified  Cox  and  Snell  bounds  will  be  compared 
to  the  unmodified  Cox  and  Snell  bounds  (the  conservative 
bound).  If  a  bound  can  be  found  that  has  adequate  coverage 
wnile  providing  tighter  (greater  rel  tightness)  upper  error 
limits  than  the  Stringer  bound,  then  that  bound  would 
provide  auditors  a  better  tool  with  which  to  audit 
accounting  populations. 

Research  Objective  3.  The  tnird  objective  is  to 
examine  the  performance  of  the  selected  bounds  by  lowering 
the  nominal  confidence  level  to  85%.  By  lowering  the 
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confidence  level,  it  may  be  possible  to  achieve  coverages 
equal  to  higher  confidence  levels  buc  resulting  in  tighter 
bounds . 

Research  Objective  4.  The  fourth  research  objective  is 
to  incrementally  loosen  each  of  the  calculated  bounds  by  5%, 
1u%,  and  20%  and  then  examine  the  resulting  bounds  for 
robustness,  coverage,  and  rel  tightness.  By  increasing  the 
bounds  in  this  manner,  it  may  be  possible  to  create  a  robust 
bound . 

Researcn  Qpjective  5.  Tne  fifth  research  objective  is 
to  examine  the  performance  characteristic  of  power  in  the 
selected  bounds  to  be  tested.  The  power  of  a  model  is  the 
percent  of  the  time,  in  repetitive  sampling,  a  bound 
correctly  predicts  that  there  is  no  material  error  in  a 
given  study  population.  The  probability  that  the  point 
estimate  from  an  upper  error  limit  model  equaling  the  true 
error  in  an  accounting  population  is  small.  An  auditor 
controls  the  risk  (Type  I  error)  that  the  model  does  not 
predict  a  material  error  in  the  population  when  there 
actually  is  a  material  error  by  the  chosen.  In  tnis  study, 
tne  performance  characteristics  of  robustness  and  coverage 
test  tne  possibility  of  a  Type  I  error.  Also  important  to 
the  auditor  and  what  this  research  objective  will  measure  is 
how  well  a  particular  oound  avoids  a  Type  II  error, 
predicting  a  material  error  in  the  population  when  tne 
actual  error  is  not  material. 
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II 


METHODOLOGY 


Introduction 

This  chapter  presents  the  sampling  approach,  upper 
error  limit  models,  the  accounting  and  study  populations, 
and  the  performance  measurement  procedures  for  the  models. 

To  remain  consistent  with  previous  research,  notation  in 
this  research  will  be  similar  to  that  used  by  Phillips 
(1985).  Much  of  the  work  in  this  research  uses  the  Stringer 
and  modified  Stringer  bounds  for  comparison  purposes.  The 
Stringer  bound  is  used  in  practice,  and  past  research  has 
shown  it  to  be  extremely  conservative.  The  main  computer 
program,  modified  for  use  in  this  research,  is  provided  in 
Appendix  A. 

The  Sampling  Approach 

Johnson,  Leitch,  and  Neter  (1981)  showed  that  dollar- 
unit  sampling  is  the  most  appropriate  sampling  method  to  use 
for  accounting  populations  [7:291J.  Dollar-unit  sampling, 
by  providing  higher  book  value  line  items  higher  probability 
of  being  selected  in  an  audit  sample,  gets  more  information 
out  of  accounting  populations  than  classical  sampling 
methods.  Therefore,  dollar-unit  sampling  will  be  used  for 
this  research.  In  dollar-unit  sampling,  the  sampling  unit 
is  the  individual  dollar.  The  following  steps  outline  the 
procedure  for  implementing  a  dollar-unit  sample.  First,  the 
book  value  of  each  line  item  in  the  accounting  population 


must  be  placed  into  a  cumulative  list  and  summed  to  provide 
the  total  reported  book  value  (T) .  A  random  sample  of  the 
proper  size  is  then  selected  from  the  dollars  one  to  T. 

Each  number  in  the  sample  refers  to  a  line  item  in  the 
population.  The  entire  line  item  to  which  the  sampled 
dollar  belongs  is  then  audited.  Lastly,  the  total  error 
found  in  an  audit  unit  is  prorated  proportionately  to  each 
dollar  in  the  line  item.  For  this  research,  a  fortran 
computer  program  will  be  modified  and  used  to  draw  500 
replications  of  systematic  random  sample  size  200  dollar- 
units  from  each  of  the  study  populations. 

Upper  Error  Limit  Bounds 

The  upper  error  limit  bound  is  a  dollar  figure  in  which 
an  auditor  would  be  confident  at  the  nominal  level  that  the 
total  dollar  error  in  a  population  would  not  exceed.  Both 
Bayesian  and  non-Bayesian  upper  error  methods  in  their 
modified  form  were  included  in  this  study.  The  non-Bayesian 
method  was  the  modified  Stringer  bound  which  was  the  basis 
by  whicn  the  modified  Cox  and  Snell  bounds,  the  Bayesian 
methods,  were  compared.  The  modified  upper  error  methods 
involve  adding  dollar  errors  identified  in  the  dollar-unit 
samples  to  the  product  of  the  bound's  upper  error  limit  and 
the  number  of  unsampled  dollars  remaining  in  the  study 
population.  This  method  eliminates  extrapolating 
probabalistic  error  limits  to  the  errors  in  the  sampled 
units  and  results  in  less  conservative  (tighter)  upper  error 
limits. 
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Tine  Stringer  Bound.  The  Stringer  method  was  introduced 
by  Kenneth  W.  Stringer  in  1963  and  uses  a  Poisson 
approximation  to  the  binomial  distribution.  The  Stringer 
upper  error  limit  can  be  broken  down  into  three  components: 
basic  precision,  most  likely  error,  and  precision  gap 
widening.  The  basic  precision  (BP)  establishes  a  minimum 
upper  error  limit  which  an  auditor  would  estimate  if  no 
errors  were  found  in  the  sample.  Each  error  found  in  a 
sample  raises  the  upper  error  limit  from  this  minimum  using 
tne  most  likely  error  and  precision  gap  widening.  The  most 
likely  error  (MLE)  is  the  auditors  best  estimate  of  the 
probable  error  in  the  population  based  on  projecting  the 
dollar-unit  sample  error  onto  the  population  [8:1 J. 

Precision  Gap  Widening  (PGW)  represents  the  amount  the  total 
precision  gap  has  widened  as  a  result  of  finding  errors  in 
the  sample.  The  sum  of  each  of  the  components  (BP+MLE+PGW) 
is  known  as  the  upper  error  limit  (UEL)  factor.  Figure  2.1 
provides  an  example  of  a  Stringer  upper  error  limit 
resulting  from  finding  four  errors  of  100%,  80%,  40%,  and 
20%  taintings  in  a  sample  of  200  dollar-units  of  a 
81,000,000  population  evaluated  at  the  90%  confidence  level. 
For  conservatism,  it  is  important  to  list  the  errors  in 
descending  order  oy  error  tainting.  This  method  gives 
larger  value  in  esimating  the  upper  error  limits  to  the 
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Basic  Precision 

2.31 

100% 

$5,000 

$1  1,550 

Most  Likely  Error 

1st  error 

1.00 

100 

5, 000 

5,  000 

2nd  error 

1 .00 

80 

5,000 

4,000 

3rd  error 

1.00 

40 

5,  000 

2,  000 

4th  error 

1  .00 

20 

5,000 

1  ,000 

Precision  Gap  Widening 

1st  error 

0.58 

100 

5,000 

2,  900 

2nd  error 

0.44 

80 

5,000 

1 ,760 

3rd  error 

0.36 

40 

5,  000 

720 

4th  error 

0.31 

20 

5,000 

310 

Total  Upper  Error 

Limit 

$29,0 

Figure  2.1  Example  of  the  Stringer  Method  at  the  90% 
Confidence  Level. 


Smith  (1979)  suggested  modifying  the  Stringer  method  in 
order  to  provide  a  tighter  bound.  The  modification 
consisted  of  applying  the  calculated  upper  error  limit 
factor  to  the  unsampled  line  items  in  the  population  and 
then  adding  the  actual  dollar  error  found  in  the  sample  to 
this  result.  The  modified  Stringer  bound  can  be  expressed 
as : 

MSTR  =  (TyUg/n)  (UEL  Factor)  +LeL  (1) 


where 


"yus  =  t'le  total  book  value  of  unsampled  line 
items 


n  =  the  sample  size 

UEL  Factor  =  the  same  UEL  Factor  as  found  in  the 


unmodified  Stringer  method 
=  the  dollar  value  of  the  error  in  the  ith 
sample  line  item  in  error. 

For  illustration  purposes,  if  the  example  in  figure  2.1 
were  such  that  the  total  book  value  of  the  sampled  line 
items  is  §100,000  and  the  book  value  of  the  line  item  with 
the  100%  taint  is  §1,000,  the  80%  taint  is  1,440,  the  40% 
taint  is  §840,  the  20%  taint  is  §480,  and  the  sample  size  is 
200  dollar-units,  the  modified  Stringer  method  could  be 
calculated  as  follows: 


MSTR  =  (§900,000/200X5.848)  +  §1,460  =  §27,776 


As  can  be  seen,  the  modification  to  the  Stringer  bound 
reduces  the  upper  error  limit  by  §1,464,  thereby  tightening 
the  bound  in  this  particular  example  by  5%. 

The  Cox  and  Snell  Bound.  Cox  and  Snell  (1979)  used 
infinite  population  theory  to  develop  a  Bayesian  upper  error 
method.  The  Cox  and  Snell  bounds  are  differentiated  by  the 
prior  parameter  selections  or  settings  of  the  dollar-unit 
error  rate  (7T)  ,  the  mean  taint  (/*t),  and  the  prior 

probability  distribution  of  7r  (designated  as  an)  and 
(designated  as  <fy)  .  The  following  calculation  is  performed 
in  developing  an  unmodified  Cox  and  Snell  upper  error  limit: 
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a  =  where  a £  is  the  expected  value  of  the 

variance  of  the  population  error  rate 

0  =  tne  given  expected  value  of  tne  mean  taint 

b  =  +  2  where  ofc  is  the  expected  value  of  the 

variance  of  the  population  mean  taint 

n  =  the  dollar-unit  sample  size 

m  =  observed  number  of  sample  dollar-units  in  error 

t  =  tne  sample  mean  taint  or  simply  the  sum  of  the 
sample  taints  divided  by  m 

FI ( 1-a)  ; 2 (m+a) , 2 (m+b) J  =  the  F  distribution  value  at 
the  nominal  confidence  level  of  (1-a)  witn  2(m+a) 
and  2 (m+b)  degrees  of  freedom  in  the  numerator  and 
denominator  respectively 

Ty  3  total  population  book  value 

For  example,  if  a  sample  of  200  dollar-units  has  four 
errors  with  taints  of  .2,  .4,  .8,  and  1.0  in  a  §1,000,000 
population,  the  resulting  Cox  and  Snell  upper  error  bound 
from  prior  parameter  settings  of  7Tq  equals  .20,  cr^  equals 

equals  .40,  ana  equals  .20,  will  be  §31,186  at  the 
90%  nominal  confidence  level.  If  another  error  were  found 
in  a  sample  of  the  same  population  with  the  same  taints  as 
the  previous  example  (meaning  a  taint  of  .6)  and  using  the 
same  prior  parameter  settings,  tne  resulting  90%  confidence 
level  bound  will  be  §35,141.  Thus,  finding  an  additional 
error  increases  the  bound  by  §1,757.  If,  in  the  first 
example,  a  different  Cox  and  Snell  bound  were  chosen  which 
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had  prior  parameter  settings  of  rrQ  equals  .20,  an  equals 
.20,  Mq  equals  .40,  and  equals  .20,  the  resulting  upper 
error  limit  bound  at  the  90%  confidence  level  would  be 
§22,280 . 

Tne  modified  Cox  and  Snell  bound  is  calculated  in  the 
same  manner  as  the  modified  Stringer  bound.  It  is  obtained 
by  statistically  projecting  the  upper  error  limit  to  the 

total  dollar  value  of  the  unsampled  line  items  (Tyus)  then 
adding  the  total  dollar  value  of  the  actual  errors  found  in 
the  sample  (Xe^)*  The  modified  Cox  and  Snell  bound  can  be 
written  as: 


MCS  =  (UELCS) (Tyus)  +  (4) 

If,  in  tne  first  example  given  for  the  unmodified  Cox 
and  Snell  bound,  a  modified  Cox  and  Snell  bound  is 
calculated  and  the  total  dollar  value  of  the  sampled  line 
items  is  §100,000,  then  the  resulting  90%  confidence  level 
bound  will  be  §29.447.  The  modification  therefore  reduced 
the  upper  error  limit  by  §1,939  or  tightened  the  bound  by 
6.2  percent. 

Prior  Probability  Parameter  Settings 

As  previously  explained,  each  Cox  and  Snell  bound  is 
distinguishable  by  the  particular  prior  probability  values 
for  the  dollar-unit  error  rate  (7 r)  ,  the  mean  taint  (^t),  and 
the  standard  deviations  of  both  parameters  (  and  ) .  The 
particular  bounds  to  be  examined  in  this  study  were  selected 


cased  upon  tne  results  of  the  Phillips  (  1985)  research  in 
which  no  modified  Cox  and  Snell  bounds  were  found  to  be 
robust.  Five  modified  Cox  and  Snell  Dounds,  eacn  somewhat 
looser  than  the  bounds  in  the  Phillips  (1985)  research,  were 
examined . 

Prior  expected  values  for  the  dollar-unit  error  rate 
(7r0)  included  in  this  research  were  four  bounds  at  .20  and 
one  at  .05  with  each  having  the  same  standard  deviation  ( (%•) 

of  .20.  Prior  expected  values  for  the  mean  taint  (/*q) 
ranged  from  small,  .05,  to  medium,  .10  and  .20,  to  large, 

.40.  The  associated  standard  deviations  (Gy)  were  selected 
to  give  variability  and  ranged  from  .10  to  .20.  Table  2.1 
provides  the  specific  prior  expectation  parameters  for  each 
of  the  modified  Cox  and  Snell  Dounds  examined  in  this 
researcn.  Note  that  the  different  bounds  are  distinguished 
by  A  through  E.  This  convention  of  identification  will  be 
maintained  throughout  this  research  and  in  every  case  refers 
to  the  particular  probability  parameter  settings  in  Table 
2.1.  In  order  to  accomplish  the  incremental  modifications 
as  specified  in  researcn  objective  4,  for  every  modified  Cox 
and  Snell  bound  calculated,  three  more  bounds  were 
calculated.  The  modified  Cox  and  Snell  bounds  were 
multiplied  by  1.05,  1.10,  and  1.20  to  provide  5,  10,  and  20 
percent  loosening. 
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Table  2.1 


Prior  Probability  Parameter  Settings 
of  Selected  Cox  and  Snell  Bounds 


Bound 

^0 

°7T 

M0 

a 

b 

AMCS 

.05 

.20 

.40 

.10 

.06 

18.00 

BMCS 

.20 

.20 

.05 

.20 

1  .00 

2.06 

CMCS 

.20 

.20 

.10 

.15 

1  .00 

2.44 

DMCS 

.20 

.20 

.20 

.20 

1  .00 

3.00 

EMCS 

.20 

.20 

.40 

.20 

1.00 

6.00 

The  Accounting  Population 

Tne  main  accounting  population  used  for  tnis  researcn 
was  acquired  from  the  University  of  Michigan  through  Dr. 
James  Godfrey.  It  consists  of  80b9  inventory  accounts  of  a 
major  pharmaceutical  manufacturer.  Account  balances  of 
$100,000  and  greater  have  been  deleted  on  the  presumption 
that  they  would  be  audited  on  a  100%  basis  in  actual 
practice.  Specific  characteristics  of  the  accounting 
population  are  listed  in  Table  2.2. 


Table  2.2 

Descriptive  Statistics  of 
Accounting  Population  Book  Values 


Number  of  Accounts 
Mean 

Standard  Deviation 

Min  imum 

Maximum 

Total 

Variance 

Skewness 

Kurtosis 

Coefficient  of  Variation 


8069 
$1,191.39 
$->,  7  .>9.53 
$0.01 
$99, 328. 32 
$9,613,301.00 
14,360,566.76 
1  1  .4292 
189.927 
3.18 


The  Study  Populations 


In  order  to  maintain  a  high  degree  of  external 
validity,  the  study  populations  were  created  with  extreme 
care.  As  noted  previously  in  the  background,  Johnson, 
Leitch,  and  Neter  (1981)  found  while  inventory  populations 
contain  both  overstatement  and  understatement  errors,  the 
majority  of  inventory  populations  still  contained 
significant  numbers  of  overstatement  errors.  Only 
overstatement  errors  were  included  in  the  population  because 
the  accounting  population  for  this  study  is  an  inventory 
account  and  Cox  and  Snell  bounds  are  not  designed  to  deal 
with  understatement  errors. 

A  broad  range  of  mean  error  rates  (7T0)were  selected  for 
this  study  to  provide  a  rigorous  test  for  the  models 
examined.  The  six  mean  error  rates  are  .80  (high),  .25  and 
.10  (medium),  and  .5,  .01,  and  .005  (low).  If  a  model 
proves  robust  (ie,  correct  approximately  (1-a)  x  100  percent 
of  the  time)  for  this  range  (favorable  to  unfavorable)  of 
mean  error  rates,  auditors  need  not  be  concerned  about 
accurately  predicting  prior  parameters  for  the  accounting 
population  they  are  auditing.  They  can  be  confident  at  the 
nominal  level  that  no  matter  what  error  rate  is  encountered 
in  the  audited  population,  the  actual  error  will  not  exceed 
the  predicted  upper  error  limit  provided  by  the  model. 

Prior  to  seeding  the  population  with  errors,  the 
population  was  sorted  from  low  book  values  to  high  book 


values.  The  population  was  then  divided  into  two  strata 
such  that  the  total  book  values  in  each  statum  were  equal 
(total  low  book  value  =  total  high  book  value).  The  low 
book  value  accounts  were  identified  as  stratum  1  and  the 
high  book  value  accounts  were  stratum  2. 

Four  different  combinations  of  mean  taints  were 
introduced  into  the  population  by  strata  as  shown  in  Table 
2.3.  They  were  selected  to  provide  significant  contrasts 
among  the  study  populations  with  a  range  from  very  high  to 
very  low  mean  taints. 

Table  2.3 

Mean  Taints 

Stratum  1  Stratum  2 

.8  .4 

.4  .2 

.2  . 1 

.1  .05 

Errors  were  seeded  into  the  population  to  reflect  the 
findings  of  Johnson,  Leitch,  and  Neter  (1981).  First,  10% 
of  all  errors  in  stratum  1  were  randomly  assigned  100% 
taints.  Second,  the  other  90%  of  the  errors  in  stratum  1 
and  all  of  the  errors  in  stratum  2  were  generated  using  an 
exponential  distribution  with  the  mean  values  shown  in  Tabl 
2.2.  In  order  to  make  the  study  populations  more  realistic 
the  exponential  taint  distributions  were  truncated  at  1.0  i 
accordance  with  Godfrey  and  Neter  (1984).  Thus, 


(6  mean 


error  rates  x  4  mean  error  taints)  different  study 
populations  were  created  to  examine  the  performance  of  the 
bounds.  The  particular  characteristics  of  each  study 
population  are  provided  in  Appendix  B. 
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Performance  Measurement  Procedures 

Performance  of  the  modified  Cox  and  Snell  bounds  was 
based  upon  robustness,  coverage,  rel  tightness  and  power 
with  respect  to  changes  in  error  rate,  mean  taint,  and  error 
amount  intensity.  A  robust  bound  must,  with  repetitive 
sampling,  be  correct  (1  -a)  x  100  percent  of  the  time  for 
all  study  populations.  Correctness  is  defined  as  a 
calculated  upper  error  limit  equaling  or  exceeding  the 
actual  error  in  the  study  population.  Coverage  is  the 
percent  of  correct  bounds  out  of  500  replications  for  each 
study  population.  The  mean  rel  tightness  measure 
constitutes  the  degree  to  which  a  calculated  modified  Cox 
and  Snell  bound  reduced  the  projected  error  in  a  study 
population  as  compared  to  that  indicated  by  the  unmodified 
Cox  and  Snell  bound.  That  is,  with  repetitive  sampling  in  a 
given  population,  mean  rel  tightness  is  the  ratio  of  the 
mean  unmodified  Cox  and  Snell  bound  to  the  mean  modified  Cox 
and  Snell  bound  and  can  be  written  algebraically  as 


replication  i. 


MCSi  =  the  modified  Cox  and  Snell  bound  for 
replication  i,  and 
i  =  1  ,  . .  . ,  500 . 

For  measurement  of  mean  rel  tightness  for  the  different 
nominal  confidence  level  bounds,  the  conservative  bound  was 
the  90  percent  nominal  confidence  level  bound  and  can  be 
written  algebraically  as 

Mean  Rel  Tightness  =  (£MCS90i/MCS85i )/500 

where 

MCS90  =  the  modified  Cox  and  Snell  bound  at  the  90 
percent  nominal  confidence  level  for 
replication  i 

MCS85  =  the  modified  Cox  and  Snell  bound  at  the  85 
percent  nominal  confidence  level  for 
replication  i 
i  -  1  ,  .  .  .  ,  500  . 

In  the  context  of  an  auditing  situation,  an  auditor 
generally  sets  a  level  of  error  in  a  population  at  or  above 
which  the  auditor  will  reject  the  population  as  being 
correct.  This  level  of  error  is  known  as  the  level  of 
materiality  or  simply  the  material  error.  When  the  true 
population  total  error  amount  is  less  than  materiality  and 
the  correct  decision  is  to  accept  the  population,  a  bound 
which  provides  adequate  coverage  may  provide  bounds  which 
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equal  or  exceed  the  material  amount  and  cause  the  auditor  to 
make  the  incorrect  decision  of  rejecting  the  population. 

The  probability  of  rejection  performance  or  power  of  a 
particular  upper  error  limit  bound  is  determined  by  the 
level  of  materiality  selected  as  criteria  for  measurement. 

In  order  to  measure  the  selected  bounds  in  this  research  for 
power,  three  levels  of  materiality  were  selected:  five 
percent,  two  and  one-half  percent,  and  one-tenth  of  one 
percent  of  the  total  book  value.  The  low  materiality  level 
of  one-tenth  of  one  percent  was  selected  for  accounting 
populations  in  which  possibly  a  single  high  value  item 
missing  might  be  of  significance  or  material.  Therefore, 
for  those  study  populations  wnich  do  not  contain  a  material 
error,  probability  of  rejection  performance  can  be  measured 
as  the  percent  of  the  time,  in  repetitive  sampling,  the 
bound  equals  or  exceeds  the  selected  level  of  materiality. 


Analysis  of  Results 


The  study  populations  were  sorted  in  ascending  order 


based  on  the  ratio  of  total  dollar  error  ( Te )  to  the  total 
population  nook  value  (Ty)<  This  ratio  is  Known  as  the 
error  amount  intensity  (EAI)  and  can  be  expressed 
algebraically  as  the  line  item  error  rate  times  the  mean 
taint  or 


EAI  =  Tli  x  Tt 

where 


(5) 


11 


v  v  T'M 


^LI  -  TyteJ/Ty  and 
Tt  *  Te/TY(e) 

Ty  =  total  book  value  of  all  line  items 

Ty(s)  =  total  book  value  of  line  items  containing  errors 
=  total  population  error  amount 

The  study  populations  were  then  divided  into  four  EAI 
categories,  low,  medium,  high,  and  very  high  based  on 
"logical  breaks"  in  the  EAI  range.  Table  2.4  identifies  the 
EAI  categories  and  their  ranges  while  Table  2.5  identifies 
tne  specific  EAI  of  each  study  population  and  the  EAI 
category  in  which  the  study  population  lies.  Each  bound 
will  be  examined  over  all  study  populations  for  its 
performance  in  robustness  and  by  each  EAI  category  for  its 
performance  in  coverage,  mean  rel  tightness,  and  power. 
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TABLE  2.4 


ERROR  AMOUNT  INTENSITY  ( EAI ) *  CATEGORIES 


Number  of  Study 

Category _ EAI  Range _ Populations 


Low 

.00022  - 

.00170 

6 

Medium 

.00442  - 

.04039 

12 

High 

.06342  - 

.12081 

4 

Very  High 

.25091  - 

.37291 

2 

Total 


Table  2.5 

Study  Populations  (P)  by  Error  Amount  Intensity 


Low 

P 

(6) 

EAI 

Me  d i um  (12) 

P  EAI 

P 

High  (4) 

EAI 

Very  High 

P  EAI 

21 

.00022 

20 

.00442 

1 

.06342 

3  .25091 

22 

.00067 

1  3 

.00457 

7 

.07626 

4  .37291 

17 

.00083 

19 

.00562 

8 

.11360 

23 

.00121 

9 

.00707 

2 

.  12081 

1  8 

.00145 

14 

.00847 

24 

.00170 

1  5 

.01364 

10  .01375 

16  .02045 

5  .02177 

11  .02608 

12  .03839 

6  .0 


Introduction 


This  chapter  contains  the  results  of  the  analysis 
performed  on  the  data  generated  from  the  research.  The 
chapter  is  presented  in  the  order  of  the  research  objectives 
and  devotes  a  major  section  of  the  discussion  to  each.  The 
first  area  discussed  is  bound  robustness  with  respect  to  all 
the  bounds  including  the  unmodified  and  modified  Stringer 
and  Cox  and  Snell  bounds.  The  discussion  tnen  transitions 
to  mean  coverage  and  mean  rel  tightness  performance  measures 
for  each  error  amount  intensity  ( EAI )  category  for  the 
unmodified  and  modified  Cox  and  Snell  bounds.  Analysis  of 
bound  performance  when  the  nominal  confidence  level  is 
lowered  from  90  percent  to  85  percent  will  follow  with  the 
discussion  of  power  presented  last. 

Due  to  the  large  number  of  bounds,  the  reader  needs  to 
understand  the  oound  nomenclature.  The  unmodified  Stringer 
bounds  are  STR90  and  STR85  meaning  the  unmodified  Stringer 
bounds  at  the  90  and  85  percent  nominal  confidence  levels, 
respectively.  The  last  two  numbers  in  the  nomenclature 
refer  to  tne  bound's  particular  nominal  confidence  level 
which  is  either  90  or  85  percent.  As  stated  in  the  previous 
chapter,  the  Cox  and  Snell  bounds  are  distinguished  by  their 
particular  prior  probability  parameter  settings.  There  are 
five  prior  parameter  settings  (A  through  E)  as  identified  in 


Table  2.1.  The  unmodified,  modified,  and  incrementally 
modified  Cox  and  Snell  bounds'  first  letter  in  the 
nomenclature  will  refer  to  these  prior  parameter  settings. 
For  example,  ACS90  refers  to  the  unmodified  Cox  and  Snell 
bound  at  the  90  percent  nominal  confidence  level  with  prior 

parameter  settings  of  7 Tq  is  .05,  <7^  is  .20,  /Jq  *40»  and  0 ^ 
is  .10.  AMCS90  refers  to  the  modified  Cox  and  Snell  bound 
at  the  same  nominal  confidence  level  and  with  the  same  prior 
parameter  settings.  AM1M9Q  refers  to  the  first  incremental 
increase  or  5  percent  loosening  of  the  modified  Cox  and 
Snell  bound  at  the  same  nominal  confidence  level  and  with 
the  same  prior  parameter  settings.  The  second  incremental 
increase  or  10  percent  loosening  of  the  bound  is  AM2M90  and 
the  third  incremental  increase  or  20  percent  loosening  of 
the  bound  is  AM3M90. 

Measurement  of  Bound  Robustness 

The  first  research  objective  was  to  examine  the  bounds 
for  robustness.  A  (1  -a)  upper  error  limit  bound  is  robust 
if  the  bound  is  correct,  with  repetitive  sampling  from  the 
population,  (1  -  ct)  x  100  percent  of  the  time.  To  meet  this 
criterion  in  this  research,  the  particular  bound  must  attain 
actual  coverage  at  the  nominal  confidence  level  or  higher 
for  all  study  populations. 

The  unmodified  and  modified  Stringer  bounds  at  both  the 
90  and  the  85  percent  nominal  confidence  levels  were  robust 
for  all  study  populations.  Table  3.1  presents  tne  actual 


mean  coverages  attained  by  each  of  the  unmodified  and 
modified  Stringer  bounds.  All  are  well  above  the  nominal 
confidence  level  which  provides  further  evidence  of  the 


conservative  nature  of  the  Stringer  bound.  No  unmodified  or 
modified  Cox  and  Snell  bounds  were  robust. 


Table  3.1 

Mean  Coverage  of  the  Stringer  Bounds 


Bound 


Mean 

Coverage 


STR90 

MSTR90 

STR85 

MSTR85 


.998 

.990 

.992 

.973 


Bound  Performance  in  Coverage  and  Mean  Rel  Tightness 

The  second  research  objective  was  to  examine  bound 
performance  in  coverage  and  mean  rel  tightness  and  to 
identify  those  bounds  which  provided  adequate  mean  coverage. 
Coverage  is  the  percent  of  the  time  in  repetitive  sampling 
the  calculated  upper  error  limit  equals  or  exceeds  the 
actual  error  in  each  study  population.  Mean  rel  tightness 
of  the  modified  Cox  and  Snell  bounds  was  measured  with 
respect  to  the  corresponding  unmodified  Cox  and  Snell  bound. 
Although  there  were  no  robust  unmodified  or  modified  Cox  and 
Snell  bounds,  some  bounds  performed  adequately  in  mean 
coverage  for  particular  error  amount  intensity  (EAI) 
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categories.  Detailed  information  on  each  bounds  performance 
in  mean  coverage  and  mean  rel  tightness  by  EAI  category  is 
provided  in  Appendix  C. 

Both  unmodified  Cox  and  Snell  bounds  D  and  E  (DCS90  and 
EC390)  performed  well  in  the  low  EAI  category.  These  two 
bounds  had  adequate  coverage  for  all  six  of  the  study 
populations  in  this  EAI  category  both  with  mean  coverages  of 
100  percent.  The  modified  Cox  and  Snell  bounds  D  and  E 
( DMCS90  and  EMCS90)  also  performed  well  in  the  low  EAI 
category.  DMCS90  and  EMCS90  adequately  covered  all  six 
study  populations  in  this  EAI  category  both  with  mean 
coverages  of  100  percent  while  providing  18.3  percent  and 
20.7  percent  (mean  rel  tightness  measures  of  118.3  percent 
and  120.7  percent)  tighter  bounds  respectively  than  the 
unmodified  Cox  and  Snell  bounds.  This  finding  indicates 
that  the  unmodified  Cox  and  Snell  bound  performs  in  a 
relatively  conservative  fashion  in  populations  with 
relatively  low  error  amount  intensities  since  the  modified 
bounds  tightened  the  unmodified  bounds  by  18.3  percent  for 
the  D  bound  and  20.7  percent  for  the  E  bound  while  not 
declining  in  actual  mean  coverage  performance. 

While  both  modified  Cox  and  Snell  bounds  D  and  E 
(DMCS90  and  EMCS90)  had  good  mean  coverage  (100  percent)  in 
tne  low  EAI  category,  DMCS90  provided  bounds  which  were 
133.1  percent  (mean  rel  tightness  measure  of  233.1  percent) 
tighter  than  EMCS90.  The  difference  in  the  two  bounds  is  in 


the  mean  taint  prior  probability  settings  (A*q).  DMCS90  nas 
a  fl o  of  .20  while  EMCS90  has  a  Uq  «40.  This  indicates 
that  the  bounds  are  very  sensitive  to  the  mean  taint  prior 
probability  parameter  settings  when  dealing  with  accounting 
populations  with  error  amount  intensities  similar  to  those 
of  the  study  populations  in  the  low  EAI  category.  Auditors, 
knowing  they  are  dealing  with  populations  with  low  EAI, 
would  likely  desire  the  tighter  bound,  DMCS90,  since  it 
provides  actual  mean  coverages  equivalent  to  EMCS90  but 
resulting  in  much  tighter  bounds. 

In  trying  to  determine  what  characteristics  of  DMCS90 
and  EMCS90  which  made  them  perform  well  in  the  low  EAI 
category,  the  prior  probability  parameter  settings  for  the 
error  rate  and  the  mean  taint  appear  to  be  significant. 
Figure  3.1  displays  the  modified  Cox  and  Snell  bound  and 
their  respective  prior  probability  parameter  settings. 

Bounds  D  and  E  are  associated  with  a  combination  of  nigh 

mean  taint,  error  rate  and  distributions  for  both  q  and  q 
relative  to  bounds  A,  B,  and  C.  While  bound  A  has  a  high 
mean  taint,  its  low  prior  probability  error  rate  setting  may 
have  detrimentally  impacted  its  mean  coverage  performance. 
None  of  the  unmodified  or  modified  Cox  and  Snell  bounds 
provided  adequate  mean  coverage  in  any  of  the  other  three 
EAI  categories  (medium,  high,  or  very  high) . 
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Figure  3.2.  Prior  Parameter  Settings  of  Selected  Cox  and 
Snell  Bounds 


Bound  Performance  at  the  85  Percent  Nominal  Confidence  Level 
The  third  research  objective  was  to  examine  the 
performance  of  the  selected  Cox  and  Snell  bounds  when  the 
nominal  confidence  level  is  reduced  from  90  percent  to  85 
percent.  As  with  the  bounds  at  the  90  percent  nominal 
confidence  level,  only  unmodified  and  modified  Cox  and  Snell 
bounds  D  and  E  provided  adequate  coverage  in  the  low  EAI 
category  and  no  bounds  provided  adequate  coverage  in  any  of 
the  three  other  EAI  categories.  Table  3.2  presents  the 
actual  mean  coverages  at  both  nominal  confidence  levels  for 
unmodified  and  modified  Cox  and  Snell  bounds  D  and  E. 
Unmodified  and  modified  Cox  and  Snell  bound  D  ( DCS85  and 
DMCS85)  provided  acceptable  mean  coverage  for  five  of  the 
six  study  populations  with  low  EAI.  Both  bounds  had  actual 
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mean  coverages  of  91.2  percent  with  the  modified  bound 
( DMCS85 )  giving  17.6  percent  (mean  rel  tightness  measure  of 
117.6  percent)  tighter  bounds  than  the  unmodified  bound 
( DCS85 ) .  The  study  population  not  adequately  covered  by 
DCS85  and  DMCS85  was  study  population  number  which  had  the 
highest  EAI  of  the  study  populations  in  the  low  EAI 
category.  Modified  Cox  and  Snell  bound  D  ( DMCS85 )  at  the  85 
percent  nominal  confidence  level  provided  actual  mean 
coverages  above  91  percent  in  the  low  EAI  category  and 
provided  23.4  percent  (mean  rel  tightness  measure  of  123.4 
percent)  tighter  bounds  than  the  same  bound  at  the  90 
percent  nominal  confidence  level. 

Unmodified  and  modified  Cox  and  Snell  bound  E  at  the  85 
percent  nominal  level  ( ECS85  and  EMC585)  adequately  covered 
all  six  study  populations  in  the  low  EAI  category  and  both 
bounds  had  actual  mean  coverages  of  100  percent.  EMCS85 
provided  a  20.5  percent  (mean  rel  tightness  measure  of 
1.205)  bound  than  the  unmodified  Cox  and  Snell  bound  (ECS85) 
at  the  same  nominal  confidence  level  and  EMCS85  also 
tightened  the  modified  Cox  and  Snell  bound  at  the  90  percent 
nominal  confidence  level  (EMCS90)  by  20.5  percent. This 
finding  is  significant  because  much  tighter  (20.5  percent) 
upper  error  limits  may  be  predicted  at  the  85  percent 
nominal  confidence  level  with  the  same  actual  confidence 
levels  as  with  the  90  percent  nominal  confidence  level. 


Table  3.2 


Actual  Mean  Coverages  of  Cox  and  Snell 
Bounds  D  and  E  in  the  Low  EAI  Category* 


Nominal 

Level 

Unmodified 

DCS 

Bounds 

ECS 

Modified 

DMCS 

Bounds 

EMCS 

90% 

100 

100 

100 

100 

8  5% 

91 

1  00 

91 

1  00 

*  all  figures  in  percent 

Although  modified  Cox  and  Snell  bounds  D  and  E  at  the 
85  percent  nominal  confidence  level  performed  well  above  the 
nominal  level  in  mean  coverage,  bound  D  provided  much 
tignter  bounds  than  bound  E.  In  the  low  EAI  category,  D 
bounds  were  138.8  percent  (mean  rel  tightness  measure  of 
238.8  percent)  tighter  than  the  E  bounds.  This  difference 
in  tightness  accounts  for  the  decrease  in  mean  tightness  by 
the  D  bound.  Therefore,  if  an  auditor  is  dealing  with 
populations  known  to  have  low  EAI  and  he  or  she  is  willing 
to  accept  slightly  more  risk,  then  the  modified  Cox  and 
Snell  bound  D  at  the  85  percent  nominal  confidence  level 
(DMCS85)  would  be  useful. 

Bound  Performance  with  the  Incremental  Modifications 

The  fourth  research  objective  was  to  incrementally 
increase  the  modified  Cox  and  Snell  bounds  by  5,  10,  and  20 
percent  in  order  to  possibly  improve  bound  performance. 
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Since  modified  Cox  and  Snell  bounds  D  and  E  performed  well 
in  the  low  EAI  category,  the  incrementally  modified  bounds  D 
and  E  also  performed  well  in  coverage.  However,  even  with 
the  third  incremental  modification  (20  percent  loosening  of 
the  modified  Cox  and  Snell  bound)  there  was  no  improvement 
in  mean  coverage  in  the  low  EAI  category.  The  mean  rel 
tightness  measure  deteriorated  steadily  for  each  incremental 
modification.  This  was  the  case  for  both  90  percent  and  85 
percent  levels  of  confidence.  Therefore,  since  the 
incremental  modifications  in  the  low  EAI  category  provided 
no  improvement  in  coverage  and  was  detrimental  to  mean  rel 
tightness,  the  incremental  modifications  to  the  modified  Cox 
and  Snell  bounds  were  not  successful  in  improving  the 
modified  Cox  and  Snell  bounds  in  the  low  EAI  category. 

The  incremental  modifications  in  the  medium  and  high 
EAI  categories  provided  no  improvement  over  the  modified  Cox 
and  Snell  bounds.  In  the  very  high  EAI  category,  the  third 
incremental  modification  to  modified  Cox  and  Snell  bound  A 
(AM3M90  and  AM3M85)  at  both  the  90  percent  and  the  85 
percenc  nominal  confidence  level  provided  good  mean 
coverage.  These  bounds  covered  both  of  the  two  study 
populations  in  the  very  high  EAI  category  at  the  nominal 
levels  or  higher  with  a  mean  coverage  well  above  the  nominal 
levels  at  94.6  percent  and  88.6  percent  respectively.  These 
bounds,  by  loosening  the  modified  Cox  and  Snell  bound  by  20 
percent,  resulted  in  a  mean  rel  tightness  measure  looser 
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tnan  the  unmodified  Cox  and  Sneii  bound  witn  the  same  prior 
parameter  settings.  The  90  percent  nominal  confidence  level 
bound  (AM3M90)  provided  a  mean  rel  tightness  measure  of  85.6 
percent  or  14.6  percent  looser  than  the  unmodified  Cox  and 
Snell  bound.  The  85  percent  nominal  confidence  level  bound 
(AM3M85)  provided  a  mean  rel  tightness  measure  of  85.3 
percent  or  14.7  percent  looser  than  the  unmodified  Cox  and 
Snell  bound  at  tne  same  nominal  confidence  level.  AM3M85 
provided  oounds  which  were  2.1  percent  (mean  rel  tightness 
measure  of  102.1  percent)  tighter  than  the  90  percent 
nominal  confidence  level  bound  (AM3M90).  Of  the  five  bounds 
studied,  bound  A  has  the  nighest  b  value  C  L/2q J  +  2)  of 
18.00  which,  combined  with  the  nighest  incremental 
modification,  may  have  been  the  two  factors  wnich  loosened 
the  bound  sufficiently  to  provide  adequate  mean  coverage. 


R 

V. 


Measurement  of  Bound  Power 

The  fifth  research  objective  was  to  measure  the 
probability  of  a  bound  rejecting  the  population  when  in 
actuality  there  was  not  a  material  error  in  the  study 
population.  The  lowest  level  of  materiality  (.1  percent  of 
total  book  value)  was  used  to  test  modified  Cox  and  Snell 
bounds  D  and  E  in  the  low  EAI  category.  At  this  level  of 
materiality,  tnree  stuay  populations  (17,  21,  and  22)  have 
actual  errors  intensities  less  than  the  material  error. 
Modified  Cox  and  Snell  bounds  D  and  E  at  botn  the  9U  percent 
and  85  percent  nominal  confidence  levels  incorrectly 
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rejected  tne  study  populations  100  percent  of  the  time. 
Further  investigation  into  the  cause  for  the  bounds'  high 
rejection  rate  revealed  that  the  possibility  of  a  lower 
limit  to  the  bounds  may  exist.  None  of  the  unmodified  or 
modified  Cox  and  Snell  bounds  studied  in  this  research  could 
provide  an  upper  error  limit  anywhere  near  the  .1  percent  of 
the  total  book  value  level  of  materiality;  all  the  bounds 
were  well  above  this  level  of  materiality.  This  result 
indicates  that  the  Cox  and  Snell  bounds  in  this  research  may 
be  too  conservative  for  populations  with  low  error  amount 
intensities  and  where  low  levels  of  materiality  are 
important . 
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IV.  Summary,  Conclusions/  and  Recommendations 
For  Future  Research 


Summary 

The  results  of  this  study  extend  research  to  find 
Bayesian  approaches  for  auditing  accounting  populations. 

Five  Cox  and  Snell  bounds  at  both  the  90  percent  and  the  85 
percent  nominal  confidence  levels  and  modifications  to 
include  known  sample  errors  in  bound  calculations  were 
examined.  Incremental  loosening  of  5,  10,  and  20  percent  of 
each  bound  was  accomplished  to  seek  bound  improvements.  The 
power  of  each  Cox  and  Snell  bound  found  to  provide  adequate 
coverage  was  also  tested.  Performance  of  the  unmodified  and 
modified  Stringer  bounds  was  also  tested  using  the  same 


study  populations. 

The  unmodified  and  modified  Stringer  bounds  were  robust 
at  botn  confidence  levels.  Two  of  the  unmodified  and 
modified  Cox  and  Snell  oounds,  D  and  E,  at  both  confidence 
levels  had  adequate  actual  mean  coverages  in  the  low  error 
amount  intensity  (EAI)  study  populations.  In  the  low  EAI 
category,  modified  Cox  and  Snell  bound  D  at  tne  90  percent 
nominal  confidence  level  provided  18. j  percent  wnile  bound  E 
provided  20.7  percent  tighter  bounds  than  tneir  correspon¬ 
ding  unmodified  bounds.  Although  bounds  D  and  E  had  similar 
actual  mean  coverages,  modified  Cox  and  Snell  oound  D,  with 
a  mean  taint  (jJq)  of  .20,  provided  bounds  which  were  13  3 
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percent  tighter  than  bound  E,  with  a  ^  Qf  .40. 

At  the  85  percent  nominal  confidence  level,  actual  mean 
coverages  of  unmodified  and  modified  Cox  and  Snell  bound  D 
exceeded  90  percent  and  reached  100  percent  with  bound  E  in 
the  low  EAI  category.  However,  significant  increases  in 
Dound  tightness  were  noted  for  both  bounds  as  the  nominal 
level  was  decreased  from  90  to  85  percent.  Modified  Cox  and 
Snell  bound  D  had  17.6  percent  tighter  bounds  and  modified 
bound  E  nad  20.5  percent  tignter  bounds  at  the  lower  nominal 
confidence  level.  Although  tne  mean  coverages  of  modified 
Cox  and  Snell  bound  D  at  the  85  percent  nominal  level  were 
lower  than  bound  E,  the  increase  in  mean  rel  tightness  of 
bound  D  was  significant  with  138.8  percent  tighter  bounds 
than  bound  E. 

The  only  bound  to  provide  adequate  mean  coverage  in 
anotner  EAI  category  was  the  third  (20  percent)  incremental 
loosening  of  modified  Cox  and  Snell  bound  A  at  both  the  90 
and  85  percent  nominal  confidence  levels  in  the  high  EAI 
category.  This  degree  of  loosening  tne  modified  Cox  and 
Snell  bound  resulted  in  bounds  which  were  more  than  14 
percent  looser  than  the  unmodified  Cox  and  Snell  bound. 

With  respect  to  power,  no  bounds  were  acceptable  at  the 
.1  percent  of  total  cook  value  level  of  materiality.  Even 
Modified  Cox  and  Snell  bound  D  with  its  increases  in  mean 
rel  tightness  failed  the  power  test  providing  100  percent 
rejection  of  populations  without  a  material  error.  Tnis 


finding  indicates  that  even  the  Cox  and  Snell  bounds 
examined  in  this  research  may  be  too  conservative  in  low  EAI 
populations . 

Conclusions 

Even  though  the  sensitivity  of  Cox  and  Snell  bounds  to 
their  prior  probability  settings  makes  prediction  of  bound 
behavior  difficult,  this  research  has  shown  that  modified 
Cox  and  Snell  bounds  D  and  E  may  provide  auditors  with  a 
better  tool  with  which  to  audit  accounting  populations  with 
low  error  amount  intensities  than  the  Stringer  bound.  In 
the  low  EAI  these  two  bounds,  D  and  E,  had  good  mean 
coverages  while  providing  significant  increases  in  rel 
tightness . 

Recommendations  for  Future  Research 

In  this  research,  mean  rel  tightness  was  measured  with 
respect  to  a  more  conservative  bound.  Further  research 
might  be  done  using  the  actual  error  in  the  study  population 
as  a  basis.  Although  this  method  of  calculating  rel 
tightness  would  only  be  applicable  in  the  laboratory 
situation  where  actual  errors  are  known,  it  may  provide  a 
broader  basis  to  compare  bounds.  Not  only  could  one  bound 
be  compared  to  another  for  tightness,  but  the  extent  of  tne 
bound's  conservatism  is  measured  directly. 

Auditors  should  consider  the  implementation  of  modified 
Cox  and  Snell  bound  D  and  E  in  accounting  populations  with 
low  error  amount  intensities. 
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Further  research  should  be  done  with  other  modified  Cox 
and  Snell  bounds  to  determine  those  that  may  provide 
auditors  with  better  upper  error  limit  bounds. 
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Appendix  A: 


************************************************************ 

C 

C  Author:  Lt  Col  J.J.  Phillips 

C 

C  Revised:  Capt  Blaine  F.  Webber 

C 

C* ********************************************************** 

^*********************************************************** 

C 

C  MAIN  PROGRAM 

C 

C*********************************************************** 

COMMON  A(  200 ) ,BMU(5) ,API(5) ,PI(5) , TMU ( 5 ) , B9 0 ( 20 0 ) 
COMMON  385(200)  ,N, SKIP, BVARY( 9000 )  ,AVARY(9  000) 

COMMON  ERRARY ( 9000 ) ,TNTARY ( 9000 ) ,TOTARY( 9000 ) 

COMMON  JZZ,CS90 ( 5 ) ,MCS90(5) ,MMC190(5) ,MMC290(5) 

COMMON  MMC390 (5) ,CS85(5) ,MCS85(5) ,MMC185(5) ,MMC285(5) 
COMMON  MMC385 ( 5 ) , IM , TOTANT, TBAR, STR90 , STR85 
COMMON  MSTR9 0 , MSTR8 5 , BBVARY (  300) ,INDEX(9000) 

COMMON  AAVARY (300) , EERARY ( 300 ) , TTNTRY ( 300 ) 

COMMON  I INDEX ( 300 ) , PGW9 0 , PGW8 5 , BP9 0 , BPS 5 , TTOTRY (300) 
COMMON  XMLE, TELI , TEDV, TBR, I GY , ISDYPP , T YSKP , TETYY 
COMMON  SIMPER, BIGER,BIGBV, I, SAMPER, SAMPBV, ITOOT 
COMMON  TE, NX, M, FN, J ,LERA( 6 ) ,TL(4) ,TH(4) , MID, TYE 
COMMON  RANARY ( 9000 ) , DSE3D ,TOT, ME , LE , CNTIT, CNTEM 
OPEN ( 1 , FILE= ' sarapl .out  ) 

CALL  FSEEK ( 1,0,2) 

OPEN( 4 , FILE= ' bound .dat1 ) 

CALL  FSEEK (4,0,2) 

OPEN(9,FILE=' info. out’ ) 

CALL  FSEEK (9, 0,2) 

REWIND  1 
REWIND  4 
REWIND  9 
WRITEO,  100) 

100  FORMAT(//1  OX,  '  BEGIN  SIMULATION . ') 

CALL  INITA 
CALL  ERRATE 
STOP 
END 

*********************************************************** 
Subroutine  Errate 

*********************************************************** 

COMMON  A( 200 ) , BMU ( 5 ) , API ( 5 ) ,PI( 5) ,TMU(5) , B90(200) 
COMMON  B85 ( 200 ) , N , SKI P , BVARY ( 9 00 0 ) ,AVARY(9000) 

COMMON  ERRARY ( 9000 ) ,TNTARY( 9000 ) , TOTARY ( 900 0 ) 
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COMMON  JZZ , CS90 ( 5 ) ,MCS90(5) ,MMC190(5) ,MMC290(5) 

COMMON  MMC 390(5) ,C385(5) ,MCS35(5) ,MMC185(5) ,MMC285(5) 
COMMON  MMC385 ( 5 ) , IM,TOTANT, TBAR, STR90 , STR85 
COMMON  MSTR90 ,MSTR8 5 , BBVARY ( 300 ) ,INDEX(9Q00) 

COMMON  AAVARY ( 300 ) , EERARY ( 300 ) , TTNTRY ( 300 ) 

COMMON  I INDEX (300), PGW9 0 , PGW8 5 , BP9 0 , BP8 5 , TTOTRY (300) 
COMMON  XMLE,TELI,TEDV,TBR, I GY, I SDYPP , TYSKP , TETYY 
COMMON  S IMPER, BIGER, BIGBV , I , SAMPER, SAMPBV, ITOOT, TYE 
COMMON  TE, NX, M, FN, J , LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 
COMMON  RANARY { 9000 ) , DSEED , TOT, ME , LE , CNTEM 
INTEGER  CNTIT, CNTEM 
REAL  LERA( 6 ) ,MID,X( 1 ) 

DOUBLE  PRECISION  DSEED 
NR=0 
1=0 

CNTIT=0 

ISDYPP=0 

CNTEM=0 

M=0 

L=0 

KK=0 

IZZ  =  0 

DSEED= 1 23457. DO 

C  These  are  the  study  population  line  item  error  rates. 

LERA ( 1 ) =  0.8 
LERA( 2 ) =  0.25 
LERA( 3 ) =  0.10 
LERA( 4 ) =  0.05 
LERA( 5 ) =  0.01 
LERA( 6 ) =  0.005 

C  These  are  the  study  population  error  taints. 

TL ( 1 ) =. 1 0 
TL( 2 ) = . 2083 
TL( 3 ) = . 8 1 3 1 
TL( 4 ) = 1 0 . 0 
TH ( 1 ) = . 0505 
TH(  2  )  =  .  1  0 
TH ( 3 ) = . 2083 
TH ( 4 )  =  . 3  1  3 1 

C  Read  the  actual  accounting  population  values  into  arrays. 
Open( unit=8,  file=,file3, , statu s= ' old ' ) 

REWIND  8 
TOTAR¥ ( 1 )=0.0 
DO  55  i=1 ,8069 
50  Read( 3 ,* ) BV,TOT, IND 
BVARY ( I ) =BV 
ERRARY ( I )=0 .0 

TOTARY ( I ) =totary ( 1-1 )+bvary(i) 

AVARY ( I ) =0 . 0 
TNTARY ( I ) =0 . 0 
INDEX( I ) = I 
55  CONTINUE 
CLOSE ( 8 ) 

I  =  1-1 


C  Calculate  the  midpoint  value  of  the  accounting  population. 
MID=Totary( I)/2 
K=1 

C  Calculate  the  number  in  the  low  book  value  stratum  (CNTIT) 
C  and  the  number  in  the  high  book  value  stratum  (CNTEM). 

57  If (MID.GE.Totary(K) ) then 
CNTIT=CNTIT+ 1 
K=  K+ 1 
GO  TO  57 
Else 

CNTEM=I-CNTIT 
END  IF 

C  Create  the  study  population  oy  combining  6  error  rates 
C  and  4  error  taints. 

Do  2000  LE=1 , 6 

C  Call  tne  subroutine  to  identify  the  line  items  with  error. 
Call  Errlin 

C  Identify  10  percent  of  the  low  book  value  items  to  receive 
C  100  percent  taints. 

IZZ  =  0 

DO  500  KM  =  1 ,  CNTIT 

IF  ( RANARY ( KM ) . EQ . 0 )  then 
GO  TO  500 
ELSE 

IZZ  =  IZZ+1 
END  IF 

500  Continue 

IZZ  =  Int(IZZ*0. 1 } 

c  Tnis  do  loop  calculates  study  populations  by  combining 
c  each  error  rate  with  4  mean  error  taints. 

DO  1000  ME=1,4 
KNT  =  0 

DO  650  KK= 1 ,1 

If ( Ranary( KK) .NE.0)then 
If ( KK.LE. CNTIT) then 
If (KNT. LE. IZZ) then 

AVARY { KK )  =  2  *  B VARY ( KK ) 

KNT=KNT+1 

Else 

550  Call  GGUBS ( DSEED, 1 , X) 

EEXP=-TL( ME) *ALOG(X( 1 ) ) 

If (EEXP.GE. 1.0)Go  to  550 
AVARY ( KK ) =BVARY ( KK ) * ( 1+EEXP) 

End  if 
Go  to  630 
Else 

560  Call  GGUBS ( DSEED, 1 , X) 

EEXP=-TH (ME) *ALOG ( X ( 1 ) ) 

If (EEXP.GE. 1 .0)Go  to  560 
AVARY (KK)=BVARY(KK) * ( 1+EEXP) 

End  if 


Else 

AVAR* ( KK ) *  BVARY ( KK ) 

End  if 

630  ERRARY ( KK ) “AVARY ( KK ) -BVARY ( KK ) 

TNTARY { KK ) =ERRARY ( KK ) /BVARY ( KK ) 

650  Continue 

ISDYPP=ISDYPP+1 
Call  Info 
CALL  STRPOF 

Print* , 1  Processing  sample  ...  1 
ERK=TH ( ME ) 

CALL  SAMPL 
Do  900  KL= 1 , I 
Avary ( KL) *0 . 0 
Errary ( KL) =0 . 0 
Tntary ( KL) *0.0 
900  Continue 

1000  Continue 
2000  Continue 
Close ( 4 ) 

END 

C* ********************************************************** 

Subroutine  Errlin 

C* ********************************************************** 

COMMON  A( 200 ) ,BMU(5) , API ( 5 ) , PI ( 5 ) ,TMU( 5 ) , B90 ( 200 ) 
COMMON  Bd5 ( 200 ) , N, SKIP, BVARY (9000) ,AVARY(9000) 

COMMON  ERRARY( 9000 ) ,TNTARY( 9000) ,TOTARY( 9000) 

COMMON  JZ Z , CS90 ( 5 ) , MCS90 ( 5 ) , MMC1 90 ( 5 ) , MMC290 { 5 ) 

COMMON  MMC39Q ( 5) , CS8 5 ( 5 ) ,MCS85(5) , MMC1 85(5) ,MMC285(5) 
COMMON  MMC 385(5) , IM,TOTANT, TBAR, STR90 , STR85 
COMMON  MSTR9 0 , MSTR8 5 , BBVARY (300) , INDEX (9000) 

COMMON  AAVARY (300) , EERARY ( 300 ) , TTNTRY (300) 

COMMON  I INDEX ( 3 0 0 ) , PGW9 0 , PGW8 5 , BP90 , BP8 5 , TTOTRY (  300) 

COMMON  XMLE, TELI , TEDV, TBR, I GY, ISDYPP, TYSKP , TETYY 

COMMON  SIMPER, BIGER,BIGBV, I, SAMPER, SAMPBV, ITOOT, TYE 

COMMON  TE, NX, M, FN, J,LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 

COMMON  RANARY ( 9000 ) , DSEED,TOT,ME,LE,CNTEM 

INTEGER  IX( 9000 ) ,G, NSAMP, NPOP, IP, MPOP, MSAMP, IER 

INTEGER  CNTIT, CNTEM 

REAL  LERA( 6 ) , POP ( 1 ) , SAMP ( 1 ) , MID 

DOUBLE  PRECISION  DSEED 

J  =  0 

K=0 

KM=0 

G=0 

NSAMP =0 
IP  =  0 
NPOP=I 
MPOP=NPOP 
DO  5  KL= 1 , I 
IX( KL) =0 
5  CONTINUE 
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Do  10  J=  1  , 1 

Ranaryi J) =0.0 
10  Continue 

C  Calculate  the  number  of  line  items  to  be  in 
C  er  ror  (  LERA(  M)  )  . 

NSAMP=I*L£RA(LE) 

C  Call  a  random  number  generator  which  gives  you  NSAMP 
C  number  of  random  numbers  in  array  IX. 

C  Transfer  the  line  items  to  be  in  error  into  array  RANARY 
20  Call  GGSRS(DSEED,0,NPOP, IP, MPOP, POP, NSAMP, MSAMP, 

*  SAMP , IX, IER) 

Do  30  IG=1 , NSAMP 
K=IX( IG) 

RANARY (K)=K 
30  Continue 
Return 
End 

C* ********************************************************** 

C 

SUBROUTINE  SAMPL 

*********************************************************** 

COMMON  A( 200 ) , BMU {  5 )  , API ( 5 ) , PI ( 5 )  , TMU ( 5 ) , B90 ( 200 ) 
COMMON  B85 ( 200 ) ,N, SKIP, SVARY( 9000 ) ,AVARY(9000) 

COMMON  ERRARY ( 9000 ) ,TNTARY( 9000 ) ,TOTARY ( 9000 ) 

COMMON  JZ  Z , CS90 ( 5 ) ,MCS90(5) ,MMC190(5) ,MMC290(5) 

COMMON  MMC390 (5) ,CS85(5) ,MCS85(5) ,MMC1 85(5) ,MMC285 ( 5 ) 
COMMON  MMC385 (5) , IM, TOTANT, TBAR, STR90 , STR85 
COMMON  MSTR90,MSTR85,BBVARY( 300 ) , INDEX (9000) 

COMMON  AAVARY (300) , EERARY( 300) , TTNTRY (300) 

COMMON  I INDEX ( 300) , PGW9Q , PGW8 5 , BP90 , BP8 5 , TTOTRY (300 ) 
COMMON  XMLE, TELI , TEDV, TBR, IGY, ISDYPP, TYSKP 
COMMON  S IMPER, BIGER, BIGBV, I , SAMPER, SAMPBV, ITOOT,TYE 
COMMON  TE, NX, M, FN, J , LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 
COMMON  RANARY ( 9000 ) , DSEED,TOT, ME, LE ,TETYY , CNTEM 
REAL  MCS90 ( 5 ) ,MMC190(5) ,MMC290(5) ,MMC390(5) ,MCS85(5) 
REAL  MMC 18  5(5) , MMC28 5 ( 5 ) , MMC385 (  5 ) 

IG  Y=0 
J=0 

DO  100  JZZ=1 ,500 
SAMPBV=BI GBV 
SAMPER=B1GER 
XMLE=0 . 0 
PGW90=0 .0 
PGW85=0 . 0 
X=0 . 0 
IM  =  0 

TOTANT=0 . 0 
ITOOT=0 

DO  10  ITO= 1 ,200 
A( ITO) =0.00 
10  CONTINUE 


b 
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TZ=TOT 
CALL  STEP1 
TOT=TZ 
IG  Y=0 

DO  20  IK= 1 , N 

SAMPER=SAMPER+EERARY( IK) 

S AMPS V=SAMPBV+BB VARY ( IK) 

IF  (TTNTRY(IK) .NE.0.0)  THEN 
ITOOT=ITOOT+1 
A( ITOOT) =TTNTRY ( IK) 

ELSE 

GO  TO  20 
END  IF 
CONTINUE 

IF  ( ITOOT. LE.1)  GO  TO  50 
CALL  SORTA(A,200,  ITOOT) 
CONTINUE 


CALL  BOUND 1 
CALL  BOUND2 
WRITE (4,60) 
WRITE (4,70) 


WRITE (4,60)  ISDYPP, JZ Z , N, TE, LERA( LE ) ,TH(ME) 

WRITE (4,70)  IM , TOT , SAMPBV , SAMPER, S IMPER 
60  FORMAT( IX, 3( 14, 2X) , 3 ( FI  2 . 2 , 2X) ) 

70  FORMAT ( 1X,I5,2X,4(F12.2, 2X) ) 

WRITE ( 4 , 80 ) STR90 , MSTR90 , STR85 , MSTR85 

WRITE( 4 , 80 ) CS90 ( 1 ) ,CS90(2) ,CS90( 3) , CS9 0 ( 4 ) , CS90 ( 5 ) 

WRITE (4,80) MCS90 ( 1 ) ,MCS90(2) ,MCS90 ( 3 ) ,MCS90 ( 4 ) , 

*  MCS90( 5 ) 

WRITE (4,80) MMC 190(1 ) , MMC1 90 ( 2 ) , MMC 1 90 ( 3 ) , 

*  MMC 1 9  0 ( 4 ) , MMC 190(5) 

WRITE (4,80) MMC290 ( 1 ) ,MMC290(2) ,MMC290(3) , 

*  MMC29  0(4), MMC2  90(5) 

WRITE (4,80) MMC390 ( 1 ) ,MMC390(2) ,MMC390(3) , 

*  MMC 390 ( 4 ) , MMC390 ( 5 ) 

WRITE (4,80) CS85 ( 1 ) ,CS85(2) ,CS85(3) ,CS85(4) ,CS85(5) 
WRITE (4,80) MCS8  5(1 ) ,MCS85(2) ,MCS85(3) ,MCS85(4) , 

*  MCS85 ( 5 ) 

WRITE (4,80) MMC 135(1 ) ,MMC1 85 ( 2 ) ,MMC1 85( 3 ) , 

*  MMC 1 8  5 ( 4 ) , MMC 185(5) 

WRITE (4,80) MMC23  5  (  1  ) ,MMC285(2) ,MMC285(  3) , 

*  MMC285 ( 4 ) , MMC285 ( 5 ) 

WRITE ( 4,80 )MMC385(1 ) ,MMC385(2) ,MMC385(3) , 

*  MMC385 ( 4 ) , MMC385 ( 5 ) 

80  FORMAT (1X,5(F14.2,2X) ) 

100  CONTINUE 
RETURN 
END 

C 

SUBROUTINE  STEP1 
C 

COMMON  A( 200 ) , BMU (5),API(5),PI(5), TMU ( 5 ) , B90 ( 200 ) 


COMMON  B85 ( 200 ) , N, SKIP, BVARY ( 9 000 ) ,AVARY(9000) 

COMMON  ERRARY ( 9000 ) /TNTARY ( 9000 ) ,TOTARY( 9000 ) 

COMMON  JZZ,CS90(5) ,MCS90(5) ,MMC190(5) ,MMC290(5) 

COMMON  MMC 3 9 0 { 5 ) ,CS85(5) ,MCS85(5) ,MMC185(5) ,MMC285(5) 
COMMON  MMC385 ( 5 ) , IM,TOTANT,TBAR, STR90,STR85 
COMMON  MSTR90 , MSTR8 5 , B BVARY ( 300 ) , INDEX(900Q) 

COMMON  AAVARY (300) , EERARY (300) , TTNTRY ( 300 ) 

COMMON  I INDEX (300), PGW9 0 , PGW8 5 , BP90 , BP8 5 , TTOTRY (300) 
COMMON  XMLE, TELI , TEDV, TBR, I GY, ISDYPP,TYSKP 
COMMON  SIMPER,  BIGER, BIGBV,  I , SAMPER, SAMPBV,  ITOOT, l'YE 
COMMON  TE, NX, M, FN, J, LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 
COMMON  RANARY ( 9000 ) , DSEED,TOT,ME,LE,TETYY,CNTEM 


14 


14 


20 

C 

C50 

C 

C 

C 

C80 

C60 


C*** 

C 


REAL  X( 1 ) 

DOUBLE  PRECISION  DSEED 
CALL  GGUBS( DSEED, 1 ,X) 

IF  (X( 1 ) .EQ.0.0)  GO  TO 
Y=X ( 1 ) *SKIP 
DO  20  IMY= 1 , N 

T I J= Y+ ( SKIP* ( FLOAT ( IMY ) -  1 .0)  ) 

CALL  SEARCH (TI J, LINEI ) 

IGY=IG Y+ 1 

BBVARY ( I GY ) =BVARY (LINEI ) 

AAVARY ( IGY ) = AVARY (LINEI) 

EERARY ( IGY) = ERRARY (LINEI ) 

TTNTRY ( IGY ) =TNTARY ( LINEI ) 

TTOTRY ( IGY ) =TOTARY (LINEI ) 

I INDEX ( IGY) « INDEX (LINEI) 

CONTINUE 

WRITE( 1 ,50) JZZ,IGY 

FORMAT( 1  OX, ' ITERATION®  ',15,'  SAMPLE  SIZE= 

DO  60  IGH= 1 , IGY 

WRITE( 1 ,80) BBVARY (IGH) , AAVARY ( IGH) , EERARY ( IGH ) , 

*  TTNTRY ( IGH ) , TTOTRY ( IGH ) ,IINDEX(IGH) 

F0RMAT(5X,5F1 2.3,15) 

CONTINUE 

RETURN 

END 

******************************************************** 


15) 


C 

C*** 


SUBROUTINE  SEARCH ( TIJ , LINE I ) 


********* 

COMMON 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


*********************************************** 

A( 200 ) , BMU ( 5 ) , API(5) , PI ( 5 ) , TMU( 5 ) , B90 ( 200 ) 

B85  (  200 ) ,N,  SKIP, BVARY (9000) ,AVARY(9000) 

ERRARY ( 9000 ) , TNTARY ( 9000 ) ,TOTARY ( 9000 ) 
JZZ,CS90(5) , MCS90 ( 5 ) , MMC 190(5) ,MMC290(5) 

MMC39G ( 5 ) ,CS8  5 ( 5 ) ,MCS85(5) ,MMC185(5) ,MMC285(5) 
MMC385 ( 5 ) , IM,TOTANT, TBAR, STR90 , STR85 
MSTR90,MSTR8 5, BBVARY ( 300) , INDEX (9000 ) 

AAVARY( 300) ,EERARY( 300) ,TTNTRY( 300) 

I INDEX (  300) , PGW9  0 , PGW8  5 , BP9  0 , BP8  5 , TTOTRY (  300 ) 
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$ 
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COMMON  XMLE,TELI,TEDV,TBR, IGY, ISDYPP, TYSKP 

COMMON  S IMPER, BIGER, BI GBV , I , SAMPER , SAMPBV , ITOOT , TYE 

COMMON  TE, NX, M, FN, J,LERA( 6 ) ,TL( 4) ,TH( 4) ,MID,CNTIT 

COMMON  RANARY ( 9000 ) , DSEED,TOT,ME,LE,TETYY , CNTEM 

INTEGER  MAXPT, MINPT, MIDPT 

MINPT®  1 

MAXPT = I 

1 0  MIDPT® { MAXPT-MI NPT ) /2+MI NPT 

IF  ( (MAXPT-MINPT) . EQ. 1 )  THEN 
LINEI=MAXPT 
RETURN 
ELSE 

IF  ( T I J . LE . TOTARY { MI DPT ) )  THEN 
MAXPT® MI DPT 

ELSE 

MINPT=MIDPT 
END  IF 
END  IF 
GO  TO  10 
END 

C 

SUBROUTINE  SORTA( A, ND, NS ) 

C 

£*********************************************************** 
REAL  A( ND) ,TEMP 
INTEGER  I , LASTS , LASTI , SSTART 
LOGICAL  INSORT 
SSTART=NS-1 
LASTS® 1 
LASTI®LASTS 
INSORT®. FALSE. 

10  CONTINUE 

IF  ( .NOT. INSORT)  THEN 
INSORT®. TRUE. 

DO  20  I=SSTART, LASTI, -1 

IF  ( A ( I ) .LT. A( 1+1 ) )  THEN 
TEMP=A( I) 

A{ I ) =A( 1+1 ) 

A( 1+1 ) =TEMP 
INSORT®. FALSE. 

LASTS® I 
END  IF 

20  CONTINUE 

LASTI=LASTS+ 1 
GO  TO  10 
END  IF 
RETURN 
END 
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C 


SUBROUTINE  STRPOP 


*********************************************************** 

COMMON  A( 200 ) , BMU ( 5 ) ,API(5) ,PI(5) ,TMU( 5 ) , B90 ( 200 ) 
COMMON  B85 ( 200 ) , N, SKIP , BVARY (9000) , AVAR* (9000) 

COMMON  ERRARY ( 9000 ) ,TNTARY( 9000 ) ,TOTARY ( 9000 ) 

COMMON  JZ 2,0390(5) ,MCS90(5) ,MMC190(5) /MMC290(5) 

COMMON  MMC390 (5) , C585 ( 5 ) ,MC385(5) ,MMC185(5) ,MMC285(5) 
COMMON  MMC385 ( 5 ) , IM,TOTANT, TBAR, STR90,STR85 
COMMON  MSTR9 0 , MSTR8 5 , MSTR7 5 , BBVARY ( 300 ) , INDEX (9  000) 
COMMON  AAVARY (300) ,EERARY( 300) , TTNTRY (300) , 

COMMON  I INDEX (300), PGW9 0 , PGtf 8  5 , BP9 0 , BP8 5TTOTRY (300) 

COMMON  XMLE, TELI , TEDV, TBR, IGY, ISDYPP,TYSKP 

COMMON  SIMPER, BIGER, BIGBV, I , SAMPER , SAMPBV, ITOOT, TYE 

COMMON  TE, NX, M, FN, J , LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 

COMMON  RANARY { 9000) , DSEED , TOT, ME , LE , TETYY, CNTEM 

N=200 

JJUMP=0 

I JUMP=0 

IBANG=0 

BIGBV=0 . 0 

I JACK= I 

BIGER=0 .0 

S IMPER=0 . 0 

TYSKP=TOT 

33  SKIP-TYSKP/N 
IBANG=0 

DO  1 0  LI P= 1 , I JACK 
IF  (  BVARY (  HP)  .GT.  SKIP)  THEN 
J=LIP 
IBANG=1 
GO  TO  34 
END  IF 

10  CONTINUE 

IF  (IBANG.EQ.0)  GO  TO  35 

34  DO  20  K= J , I JACK 

I JUMP=I JUMP+ 1 
JJUMP=JJUMP+1 
BIGBV=BIGBV+BVARY ( K) 

BIGER=BIGER+ERRARY( K) 

S IMPER=S IMP ER+ ERRARY ( K ) 

20  CONTINUE 

TYSKP=TOT- BIGBV 
N=N- J JUMP 
I JACK3 I-I JUMP 
JJUMP=0 

IF  ( IBANG.GT. 0 )  GO  TO  33 

35  RETURN 
END 
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no  on  on 


C 

SUBROUTINE  INFO 


COMMON  A(200) , BMU ( 5 ) , API ( 5 ) ,PI(5) , TMU (5) , B9Q ( 200 ) 
COMMON  B85 { 200 ) , N, SKIP, BVARY (9000) ,AVARY(9000) 

COMMON  ERRARY ( 9000 ) ,TNTARY( 9000 ) ,TOTARY( 9000 ) 

COMMON  JZ Z , CS90 ( 5 ) ,MCS90(5) ,MMC190(5) ,MMC2  90(5) 

COMMON  MMC390 (5) ,CS85(5) ,MCS85(5) , MMC1 85(5) , MMC285 ( 5 ) 
COMMON  MMC385 ( 5 ) , IM, TOT ANT, TBAR, STR90 , STR85 
COMMON  MSTR9 0 , MSTR85 , BBVARY ( 300 )  ,INDEX(9000) 

COMMON  AAVARY( 300) , EERARY ( 300 ) , TTNTRY( 300) 

COMMON  I INDEX { 30 0 ) , PGW9 0 , PGW8 5, BP90, BP8 5 , TTOTRY (300) 

COMMON  XMLE, TELI , TEDV,TBR, IGY, ISDYPP,TYSKP 

COMMON  SIMPER, BIGER,BIGBV, I, SAMPER, SAMPBV, ITOOT,TYE 

COMMON  TE, NX, M, FN, J, LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 

COMMON  RANARY ( 9000 ) , DSEED , TOT, ME, LE , TETYY , CNTEM 

NX=0 

TE=0 . 0 

TOT=0 . 0 

M=0 

TYE=0 . 0 
TELI  =  0 . 0 


TEDV-0 . 0 
TBR=0 .0 
TETYY=0 . 0 
DO  1  JACK* 1,1 
NX=NX+1 

TE=TE+ERRARY ( JACK ) 

TOT-TOT+BVARY ( JACK ) 

IF  (ERRARYC JACK) .NE. 0.0)  THEN 
M=M+1 

TYE=T YE+BVARY ( JACK ) 

END  IF 

1  CONTINUE 

TELI-FLOAT( M) /FLOAT( NX) 

TEDV=T YE/TOT 

TBR=TE/TYE 

TETYY=TE/TOT 

WRITE (9,20)  ISDYPP, M, NX , TOT, TYE 
WRITE (9,21 )  TE, TELI, TEDV,TBR, TETYY 

20  FORMAT( 1X,3(I5,2X),2(F14.4,2X)) 

21  FORMAT ( 1X,F14.2,4(F14.6,1X) ) 

RETURN 


I 

i 


« 

i 

B8 5(200) , N, SKIP, BVARY (9000) , AVARY ( 9000 )  j 

ERRARY ( 9000 ) ,TNTARY( 9000 ) ,TOTARY ( 9000 )  i 

JZ  Z , CS90 ( 5 ) , MCS90 ( 5 ) ,MMC 190(5) ,MMC290(5)  ! 

MMC390 ( 5 ) ,CS85 ( 5 ) ,MC535(5) ,MMC185(5) ,MMC285(5) 

MMC385 ( 5 ) , IM,TOTANT, TBAR, STR90 , STR85 
MSTR9 0 , MSTR8  5 , BBVARY (  300), INDEX (9000) 

AAVARY (300) , EERARY ( 300 ) , TTNTRY ( 300 )  1 

I INDEX( 300 ) , PGW9  0 , PGW8  5 , BP9  0 , BP8  5 , TTOTRY (  300 ) 

XMLE, TELI , TEDV, TBR, I GY, ISDYPP,TYSKP 
S IMPER, BIGER, BIGBV , I , SAMPER, SAMPBV , ITOOT , TYE 
TE , NX, M, FN, J , L£RA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 
RANARY ( 9000 ) , DSEED ,TOT, ME , LE , TETYY, CNTEM 
REAL  STR90,  STR85 ,  MSTR90 ,  MSTR85 
BP90=2. 31 
BP85= 1 .90 
DO  10  IZ=1,N 

XMLE=XMLE+A( IZ) 

PGW90=PGW90+(B90( IZ)*A( IZ) ) 

PGW85=PGW85+(B85(IZ)*A(IZ) ) 

10  CONTINUE 

FN=FLOAT(N) 

STR90=( BP90+XMLE+PGW90 ) *TYSKP/FN+S IMPER 
STR85=( BP35+XMLE+PGN85 ) *TYSKP/FN+S IMPER 
MSTR90= ( TOT-SAMPBV } * ( BP90+XMLE+PGW9Q ) /FN+SAMPER 
M5TR85= (TOT-SAMPBV) * ( BP85+XMLE+PGW85 ) /FN+SAMPER 
RETURN 
END 

c*********************************************************** 

C 

SUBROUTINE  BOUND2 
C 

0*********************************************************** 

COMMON  A( 200 ) , BMU ( 5 ) , API ( 5 ) ,PI(5) , TMU ( 5 ) ,B90(200) 

COMMON  B85 ( 200 ) ,N, SKIP, BVARY( 9000 ) ,AVARY(9000) 

COMMON  ERRARY( 9000 ) ,TNTARY( 9000) ,TOTARY( 9000)  . 

COMMON  JZZ,CS90(5) ,MCS90(5) , MMC1 90(5) , MMC290 ( 5 )  j 

COMMON  MMC390 ( 5 ) ,CSd 5 ( 5 ) , MCS8 5 ( 5 ) ,MMC185(5) ,MMC285(5) 

COMMON  MMC385 ( 5 ) , IM , TOTANT, TBAR, STR90 , STR85 
COMMON  MSTR90 , MSTR8  5 , BBVARY ( 300 ) ,INDEX(9000) 

COMMON  AAVARY (3 00) , EERARY (3 00) , TTNTRY (300) 

COMMON  I INDEX (  300 ) , PGW9 G,PGW85,BP90,BP85 , TTOTRY (  3  00  ) 

COMMON  XMLE, TELI, TEDV, TBR, IGY, ISDYPP,TYSKP  i 

COMMON  SIMPER, BIGER, BIGBV, I , SAMPER , SAMPBV , ITOOT , TYE 
COMMON  TE, NX, M, FN , J , LERA( 6 ) ,TL(4) ,TH(4) ,MID,CNTIT 
COMMON  RANARY ( 900 0 ) , DS EED , TOT, ME, LE , TETYY, CNTEM 
REAL  P,D1 , D2 , X, FIRST ( 5 ) ,MCS90(5) ,MMC190(5) ,MMC290(5) 

REAL  MMC390 ( 5 ) ,MCS85(5) ,MMC185(5) ,MMC285(5) ,MMC385(5) 

INTEGER  IER 
FN=FLOAT( N) 

DO  12  IW= 1 , N 

IF  ( A( IW) .GT.0.0 )  IM= IM+ 1 
TOTANT=TOTANT+A( IW) 
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CONTINUE 

IF  (IM.EQ.O)  THEN 
TBAR=0 . 0 
ELSE 

TBAR=TOTANT/ FLOAT ( N ) 

END  IF 

DO  13  IG= 1 , 5 

FIRST{ IG) = ( { ( { FLOAT( IM ) *TBAR) + 

(BMU(IG)-I )*TMU(IG) )/ 

( (API(IG)/PI(IG) ) +  FN ) ) * ( FLOAT ( IM ) + AP I ( I G ) )/ 
( FLOAT ( IM) +BMU{ IG) ) ) 

D 1 =2* ( FLOAT ( IM) +API ( IG ) ) 

D2=2*( FLOAT ( IM)+BMU( IG) ) 

P= .  90 

CALL  MDFI  ( P , D 1 , D2 , X, IER) 

CS90 { IG) =FIRST( IG) *X*TYSKP+S IMPER 
:S90 ( IG) *FIRST( IG) *X* ( TOT-SAMPBV ) +SAMPER 
MMC1 90 ( IG) =MCS9  0 ( IG) *  1 . 05 
MMC290 ( IG) sMCS90 ( IG) *  1 .10 
MMC390 ( IG) =MCS9Q ( IG) *  1 .20 
P= .  85 

CALL  MDFI  {  P ,  D 1 , D2 , X, I ER) 

CS85(IG)=FIRST(IG) *X*TYSKP+S IMPER 
MCS85 ( IG ) =FIRST( IG) *X* ( TOT-SAMPBV ) +SAMPER 
MMC 1 8  5 ( IG ) =MCS8  5(IG)*1.05 
MMC285 ( IG) sMCS8  5 ( IG) *  1 . 10 
MMC38  5 ( IG) =MCS85 ( IG) *  1 .20 
CONTINUE 
RETURN 


SUBROUTINE  INITA 


COMMON  A(  200 ) ,BMU(5) ,API(5) ,PI(5) , TMU ( 5 ) , B90 ( 20 0 ) 
COMMON  B85 ( 200 ) ,N, SKIP , BVARY ( 9 000 ) ,AVARY(9000) 

COMMON  ERRARY (9000) , TNTARY ( 9 00 0 ) ,TOTARY ( 9000 ) 

COMMON  JZZ,CS90(5) ,MCS90(5) , MMC 190(5) ,MMC290(5) 

COMMON  MMC390 (5) , CS85 ( 5 ) , MCS85 ( 5 ) ,MMC185(5) ,MMC285(5) 
COMMON  MMC 3 85(5) , IM , TOTANT , TBAR, STR90 , STR85 
COMMON  MSTR9Q , MSTR85 , B BVARY ( 300 ) , INDEX (9000) 

COMMON  AAVARY ( 300) , EERARY ( 300 ) ,TTNTRY( 300) 

COMMON  I INDEX ( 300) , PGW9 0 , PGW85 , BP90 , BP8 5 , TTOTRY ( 3 00 ) 

COMMON  XMLE, TELI , TEDV, TBR, IGY, ISDYPP,TYSKP 

COMMON  SIMPER, BIGER, BIGBV, I , SAMPER, SAMPBV, ITOOT, TYE 

COMMON  TE, NX, M, FN, J , LERA( 6 ) ,TL(4) ,TH{4) ,MID,CNTIT 

COMMON  RANARY( 9000 ) , DS EED ,TOT, ME , LE , TETYY , CNTEM 

B  VALUES 

BMU ( 1 )-1 8.00 

BMU( 2 ) =2 .0625 

BMU( 3 ) =2 . 4444 


V  JV  _Vv  V 


BM  U  (  4  )  =  3 . 0  0 
BMU ( 5 ) =6 . 00 
A  VALUES 
API( 1 )= . 0625 
API { 2 ) = 1 .00 
API ( 3 ) = 1 .00 
API(4) =1 .00 
API ( 5 ) = 1 .00 
ERROR  RATES 
PI ( 1 )=.05 
PI( 2) =.20 
PI  ( 3 )  =  . 20 
PI ( 4 ) = . 20 
PI  ( 5 )  =  . 20 
MEAN  TAINTS 
TMU { 1  )  =  .40 
TMU ( 2 ) = .05 
TMU ( 3 ) = . 10 
TMU ( 4 )  =  . 2  0 
TMU ( 5 ) = . 40 

PRECISION  GAP  WIDENING  VALUES:  90%  CONFIDENCE  LEVEL 

B90 ( 1 )  = .  5  8 

B90 ( 2 ) = . 44 

B90 ( 3 ) = . 36 

B90 ( 4 ) = . 3 1 

B90 ( 5  )  = . 28 

B90 ( 6 ) = . 26 

B90 ( 7 ) = . 

B90 ( 8 ) = . 22 
B90 ( 9  )  =  . 2 1 
B90 ( 1 0 ) = . 20 
B90( 1 1 )=. 1 9 
B90 ( 12)=. 19 
B90 ( 1 3 ) * . 1 7 
B90 ( 1 4 )  =  .  1  7 
B90( 15)=. 17 
B90 ( 1 6 ) = . 1 6 
B90( 1 7 ) =. 1 5 
B90( 18)  =  .  15 
B90 { 1 9 ) = . 1 5 
B90 ( 20 ) = . 1 4 
B90 ( 2 1 ) = . 1 4 
B90 ( 22 ) = . 13 
B90( 23 ) =. 1 4 
B90 ( )  =  .  1  3 
B90 ( 25 ) = . 1 3 
B90 ( 2  6 )  =  .  1  4 
B90 ( 27 ) = . 1 2 
B90 ( 28 ) = . 12 
B90( 29 ) =. 1 2 
B90 ( 3  0 )  =  . 12 
B90 ( 3 1 ) = . 1 2 


SA ^ 


Ai. 
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B90 ( 32  ) 3 . 1  2 
B90 ( 3  3  )  = .  1  2 
B90 ( 34 ) = . 12 
B90 ( 35 ) =. 1 2 
B90 ( 36 ) 3 . 1 1 
B90 ( 37 ) 3 . 1 1 
B90 ( 3  8 ) 3 . 1  1 
B90 ( 39 ) 3 . 1  1 
DO  210  JC 1 =40 , 50 
B90{ JC1 )=. 10 
10  CONTINUE 

DO  220  JC2=5 1 ,59 
B90 ( JC2 ) = . 09 
20  CONTINUE 

DO  230  JC3=60 , 7 9 
B90 ( JC3 ) = . 08 
30  CONTINUE 

DO  0  JC4=30 , 94 
B90 ( JC4 ) = . 07 
CONTINUE 

DO  250  JC5=95 , 1 10 
B90 ( JC5 )  =  .0  6 
50  CONTINUE 

DO  260  JC6= 1 11,125 
B9  0 ( JC6 )  = . 0  5 
50  CONTINUE 

DO  270  JC7=1 26,140 
B90 ( JC7 ) =. 04 
70  CONTINUE 

DO  280  JC8=1 4 1 ,155 
B9  0 ( JC8 ) = . 0  3 
30  CONTINUE 

DO  290  JC9=1 56 , 170 
B90 ( JC9 ) =.02 
90  CONTINUE 

DO  300  JC1 0=171,185 
B90( JC10)=.01 
00  CONTINUE 

DO  310  JC 1 1=186,200 
B90 ( JC1 1 ) =0 . 0  0 
10  CONTINUE 

PRECISION  GAP  WIDENING  VALUES: 
B85( 1 ) =.48 
B85 ( 2 ) = . 35 
B85 ( 3 )  = . 2  9 
B85  (  4 )  =  . 25 
B85 { 5  )  =  .  23 
B85 ( 6 ) = . 21 
B85 ( 7 ) =. 1 9 
B85 ( 8 ) = . 1 8 
B85 ( 9 ) = . 1 7 
B85 ( 10)=. 17 


85%  CONFIDENCE  LEVEL 


B85 ( 1  1  )  =  .  15 
B8 5(11)=. 15 
B85( 12)  =  . 14 
B85( 13)  =  . 14 
B85 ( 1 4 )  =  .  1  3 
B85( 15)=. 1 3 
B85( 16)  =  .  13 
B85 ( 1 7 ) = . 1 2 
B85( 18)  =  .  12 
B8 5(19)  =  . 11 
B85 ( 20 )  = .  1  1 
B8  5 ( 21 )  =  .  1  1 
B85 ( 22 ) = . 1 1 
B8  5 ( 2  3 )  = .  1  1 
B85 ( )  =  .  1  1 
B85 ( 25 ) =. 1 0 
B85 ( 26 )  = .  1  0 
B85 ( 27 ) = . 1 0 
B85 ( 28 )  =  .  1  0 
B8 5(29)  =  . 10 
DO  110  JA1 =30 , 40 
B85( JA1 )  = . 0  9 
CONTINUE 
DO  120  JA2=4 1 ,49 
B85 ( JA2 ) = .08 
CONTINUE 
DO  130  JA3=50 ,69 
B85 ( JA3 )  =  . 0  7 
CONTINUE 
DO  140  JA4=70 ,89 
B85 ( JA4 ) = .06 
CONTINUE 

DO  150  JA5=90 ,105 
B85 ( JA5 )  =  . 0  5 
CONTINUE 

DO  160  JA6= 1 06 , 125 
B85 ( JA6 )  = .  0  4 
CONTINUE 

DO  170  JA7=1 26,140 
B8  5 ( JA~ )  =  . 0  3 
CONTINUE 

DO  180  JA8=1 41,155 
B8  5 ( JA8 )  = . 0  2 
CONTINUE 

DO  190  JA9=1 56 , 175 
B85(JA9)=.01 
CONTINUE 

DO  200  JA1 0=176, 200 
B85 ( JA1 0 ) =0 . 00 
CONTINUE 
CONTINUE 
RETURN 
END 
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Appandix  B:  Study  Population  Characteristics 
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.16 

814,558 
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14 
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20 

.01 
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.06 

34,431 

2,086 
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